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ABSTRACT 

We characterize the radial and angular variance of the Hubble flow in the COMPOSITE sample 
of 4534 galaxy distances. Independent of any cosmological assumptions other than the existence of 
a suitably averaged linear Hubble law, we flnd with decisive Bayesian evidence {hiB ^ 5) that the 
Hubble constant averaged in spherical radial shells is closer to its global value when referred to the rest 
frame of the Local Group rather than to the standard rest frame of the Cosmic Microwave Background 
(CMB) radiation. Angular averages reveal a dipole structure in the Hubble flow variance, correlated 
with structures within a sphere of radius 30 - 62 /i~^Mpc. Furthermore, the angular map of 
Hubble flow variance is found to coincide with the angular map of the residual CMB temperature dipole 
in the Local Group rest frame, with correlation coefficient —0.92. This suggests a new mechanism 
for the origin of the CMB dipole: in addition to a local boost it is generated by differences in the 
distance to the surface of last scattering, of a maximum ±0.35 /i~^Mpc, which arise from foreground 
structures within 65/i~^Mpc, a 0.5% effect. The dipole feature is accounted for by our position 
in a filamentary sheet between nearby voids and walls, producing a foreground density gradient on 
scales up to 65/i~^Mpc on opposite sides of the sky. This result potentially eliminates problems 
of interpretation of "bulk fiows". Furthermore, anomalies associated with large angles in the CMB 
anisotropy spectrum, and also the dark fiow inferred from the kinematic Sunyaev-Zel'dovich effect on 
small angular scales, need to be critically re-examined. 

Subject headings: cosmology: observations — cosmology: theory — distance scale — large-scale struc- 
ture of universe — cosmic background radiation — galaxies: kinematics and dy- 
namics 



1. INTRODUCTION 

It is usually assumed that the cosmic mi c rowave back- 
groun d (CMB) dipole (jKogut et all 119931 : iFixsen et aH 
Il996f ) is generated entirely by our own local peculiar ve- 
locity. A local boost by the opposite velocity then defines 
the cosmic rest frame in which we can be considered to 
be comoving observers in the background geometry of 
homogeneous isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) model. Indeed measurements of cosmo- 
logical redshifts are routinely transformed to the CMB 
frame. According to the implicit assumptions made in 
such a transformation, the cosmic rest frame so defined 
should also be the frame in which the Hubble fiow is most 
uniform, with minimal statistical variations as compared 
to other choices of the standard of rest. 

Our understanding of the Hubble fiow is, however, 
greatly comphcated by the fact that the universe is 
not completely homogeneous. Rather it appears to 
only be homogeneous in some statistical sense, when 
one averages on scales >100/i~^Mpc, the transition 
scale still being a matt er of debate (|Hogg et al.l 120051 : 
iSvlos Labini et al]l2009f ). At scales below or comparable 
to the scale of statistical homogeneity a complex pattern 
of variance in the Hubble fiow is observed. In the stan- 
dard manner of thinking about the problem, Hubble fiow 
variance is interpreted as a field of peculiar velocities of 
galaxies with respect to the expansion law of a FLRW 
model, which is linear on scales up to 2; ~ 0.1, well above 
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iWevant et al.l (pMTl ): 

of peculia r velocities 



the scale of statistical homogeneity. The CMB rest frame 
sets the standard of rest for a comoving observer at 
our location in defining the leading order linear Hubble 
law. A great deal of observational effort has gone into 
understand ing the nearby peculi ar motions so derived; 
see, e.g., 'Hud son et al.l (l2004i) : iKocevski fc EbelingI 
(I2006D: IS pringob et all () 2MZl ); iTullv et al.l (120081) : 
[Karachentsev et al. f2009|}; iLavaux et al.l (l2010li: 
Antoniou & PerivolaroDoul(3 (120101); IColin et al.l 
(2011); Dai et al. (201l]) 
Iwata fc Chamaraux ( 201 X). 

A number of studies 

(jWatkins et all 120091 : |Feldman et^ l2010f ) have found 
results which indicate persistent bulk ffows extending to 
very large scales, and which are potentially in conffict 
with the expectations of the perturbed FLRW model 
that underlies the standard Lambda Cold Dark Matter 
(ACDM) cosmology. Different data sets and methods 
of analysis produce differ ent, sometimes contradictory, 
results. For example, iTurnbull et al.l (|2012f ) made 
a recent study using 245 type la supernova (Snela) 
distances. They find a bulk fiow 249 ± 76kms~^ in 
the direction I = 319° ± 18°, 6 = 7° ± 14°, which is 
consistent with the predictions of the ACDM model, 
but also marginally consistent with the larger bulk fiow 
of 407 ± 81kms-i toward i = 287° ± 9°, 6 = 8° ± 6° 
found in the large COMPOSITE data set of mainly 
non- Snela galaxy distances studied b y I Wat kins et al.l 
(Hop). The results of Turnbull et al.l ()2012D appear to 
be inconsis tent, howe v er, with the la rger bulk fiow found 
by Kashli nskv et all (j2009al . I2010f ) using the kinetic 
Sunyaev-Zel'dovich effect. 
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In the above papers and el sewhere in the litera- 
ture, with very few exceptions (jMcClure fc Dverll2QQ7l : 
iLi fc Schwar"a[2QQ8l) . variance in the Hubble flow is gen- 
erally attributed to peculiar velocities whose radial com- 
ponents are deflned as deviations from a linear Hubble 
law 

^pec = - H^r (1) 

where z is the redshift, c the speed of light and r an ap- 
propriate distance measure. Such a deflnition implicitly 
makes a strong assumption about spacetime geometry, 
namely on the scales of interest spatial curvature can be 
neglected and the redshift associated with the Hubble 
expansion can be treated in the manner of a recession 
velocity as in special relativity. 

From the point of view of general relativity, without 
any a priori cosmological assumptions, such an assump- 
tion must be questioned. In general relativity in an arbi- 
trary spacetime background the only velocities that are 
uniquely related to observables are those corresponding 
to local boosts at a point. Given that the dust approxi- 
mation is not rigorously deflned for the complex cosmic 
web of voids and walls that constitute the prese nt day 
universe on < lOO/^-^Mpc scales (I Wiltshirel[2QTTl ). then 
the extension of the concept of a velocity in ([1]) over the 
vast distances over which space is expanding is merely 
an ansatz, whose validity remains to be justifled. In any 
general relativistic framework there must be some local 
peculiar velocities, which arise from the local dynam- 
ics of galaxies within bound clusters. However, there 
is no a priori reason for assuming that all redshifts on 
scales < 100 h~^Mpc can be treated in terms of a simple 
Doppler shift in Euclidean space, which is in practice the 
method of analysis adopted by most observationalists. 

From the point of view of general relativity, variance 
in the Hubble flow is more naturally viewed as the dif- 
ferential expansion of regions of different local densities 
which have experienced different local expansion histo- 
ries over the billions of years that have elapsed since the 
density fleld was close to uniform. While particular geo- 
metrical assumptions would lead to the standard picture 
of bulk flows in a FLRW geometry with Newtonian per- 
turbations, the lack of convergence of bulk flows to the 
CMB dipole and the puzzle of several possible anomalies 
associated both with bulk flows and the large angle mul- 
tipoles of the CMB anisotropy spectrum (Tegmark et al] 
[20031: ILand &: M agueiio 120051: Eriksen et al. 2007; 
iKashlinskv et al.l i2009a: Kim Nas elskv 2010). suggest 
that one should reconsider the problem from flrst princi- 
ples. 

In this paper we will therefore reanalyse the 
largest availab l e dat a set, th e COMPOS ITE sample of 
IWatkins et all (I2009D and Feld man et all f2010), from a 
fresh perspective. While the particular analysis we adopt 
is one which is nat urally su ggested by the cosmological 
model of lWiltshird (j2007al lbl. l2008l 1200911 , our actual anal- 
ysis is independent of any cosmological model assump- 
tions other than the most elementary one that a suit- 
ably deflned average linear Hubble law exists. We will 
flrst consider radial averages, and then angular averages. 
We then suggest a new mechanism for generation of the 
CMB dipole, consistent with our results, and discuss its 
implications. 



2. RADIAL AVERAGES 

We adopt the point of view that on scales of order 
10/i~^Mpc - 30/i~^Mpc the regional expansion history, 
and the regional average Hubble law should be deter- 
mined primarily by the relevant regional average density. 
From the point of view of any observer, under dense voids 
will appear to be expanding faster than denser wall re- 
gions, on account of the wall regions having decelerated 
more. This is true independently of whether or not there 
is a homogeneous dark energy which acts to accelerate 
the expansion by the same amount in all regions. 

The largest typical voids are shown by surveys t o have 
a diameter -30/i-^Mpc (Hoyle & Vogejeyl |2002|, 120041 : 
Pan et al. 2 011f ). and our galaxy is located in a fllamen- 
tary sheet on the edge of a Local Void of at least this 
diameter, formed from a complex of three smaller voids 
(jTullv et al.]l2008[ ). 

Although the expansion rate just to the other side of 
a local void (wall) will appear faster (slower) than av- 
erage, whole sky spherical averages that include many 
structures in different directions can be expected to have 
a reduced variance as compared to measurements in par- 
ticular directions. Furthermore, once one also averages 
on radial scales a few times larger than the largest typi- 
cal structures then the variance of the Hubble parameter 
averaged in spherical shells will reduce to a level con- 
sistent with individual measurement uncertainties. This 
provides an operational deflnition of the scale of statisti- 
cal homogeneity independent of any detailed cosmological 
model assumptions. 

A study of the spherically averaged Hubble flow, 
as a function of radial distance was undertaken by 
ILi &: Schwar3 (|2008 ) (henceforth LS08), using a sub- 
set of 54 distances fr om the HST Key project data 
(jFreedman et al.l l2001f ). Fi gure 2 of LS08 shows the 
radial variation SH{r) = {H{r) — Hq)/Hq that results 
from such an analysis, with data restricted to the range 
^min < r < Tmax, whcrc rmin = 22.5 /i~^Mpc and 
^max = 130/i~^Mpc. Furthermore, their flt is com- 
puted for redshifts referred to a single linear Hubble law 
cz = H{r)r within a sphere of radius r, as r is varied, 
and with redshifts referred to the CMB frame. 

Our flrst aim here is to perform a similar analy- 
sis to LS08 using the COMPOSITE sample of cluster 
group and galaxy distances compiled by Watkins et al.l 
(I2009D fhencef o rth W FH09) and slightly updated by 
iFeldman et all (|2010^ (henceforth FWHIO). For each 
sample object, redshift, galactic latitude and longitude, 
distance, and distance uncertainty are given. Distance 
uncertainties are about 15% for most individual galaxies. 
We include all 4,534 data points outside the local group 
r > 2 /i~^Mpc, which extends the range of distances con- 
sidered in LS08 to both smaller and larger values. 

The data in the COMPOSITE sample combines nine 
independent, full-sky samples, nearly all major peculiar 
velocity surveys published to date. Although each sur- 
vey uses a different distance measurement methodology, 
all of the surveys w ere shown to be stat istically consis- 
tent with each othe r (IWatk ins et al.ll2"009f ) . The survey of 
iLauer fc PostmanI (1994) wa s not inclu ded, as it gave in- 
consistent results (Watkins et al.l l200'9l ) . WFH09 provide 
a detailed discussion of the issues involved in combining 
subsamples with different characteristic depths and sky 
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coverages. For our analysis it is important to note that 
outside the Zone of Avoidance (ZoA) the COMPOSITE 
sample has good all sky coverage, as is seen from Fig. 1 
of FWHIO, and in Figure [T] below. 



propagation for the linear fit (|2]) in the 5th shell 
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2.1. Methodology 

Our analysis will feature two key differences to that of 
LS08. Firstly, rather than simply performing the analy- 
sis in the CMB frame, we also perform the analysis in the 
local group (LG) and local sheet (LS) frames. A com- 
parison of these frames is motivated by the Cosmological 
Equivalence Principle (|Wiltshirel l2QQ8f ): the LG frame 
corresponds to our "finite infinity region" (jWiltshird 
l2QQ7a[ ). and should be close to the frame in which the 
variance in the Hubble flow is minimized i n the ap- 
proach to cosmological averages advocated by ' Wiltshir^ 
(|2QQ7al E [2008 ). The LG has a smal l peculiar velocity 
of 66 ± 24kms-i relative to the LS (jTullv et al.ll2008f ) 
within which a local Hubble flow is first defined. 

Secon dly, on acco u nt of the small number of data 
points ILi fc Schwarzl ("2008) included all data within a 
radius r, as r was varied in steps between rmin and 
^max- This has the effect that each binned data point 
shown in their Figure 2 is correlated with the previ- 
ous data point. With 4,534 available data points in the 
COMPOSITE sample such correlations can be avoided 
by the following technique: we will minimize the sum 

Xs = [a~^{ri —cZi/H)] with respect to as a 
means of fitting a Hubble law by a standard linear re- 
gression (Press et al. 1986, Sec. 14.2), in successive in- 
dependent radial shells < r < rg+i. We consider the 
linear Hubble law with r as a function of z since all un- 
certainties have been included as distance uncertainties 
in the COMPOSITE sample. The value of the Hubble 
constant Hs computed for the sth shell is then 



and (ii) a zero point uncertainty 



s * ^ ' 

' s 



(4) 



H, 




(2) 



where cr^ denote individual distance uncertainties in 
/i~^Mpc. 

The total uncertainty for Hs in each shel0, , is deter- 
mined by adding the following uncertainties in quadra- 
ture: (i) the uncertainty determined by standard error 



^ While the measurement uncertainties in redshifts are neghgible, 
using the standard pecuhar velocity framework a uniform velocity 
noise uncertainty was added in quadrature to H^ai in FWHIO in 
defining the maximum likelihood weights. In the peculiar velocity 
framework galaxy motions are modeled using linear theory. The ve- 
locity noise term then accounts for the fact that individual peculiar 
velocities deviate from the local value of the linear peculiar velocity 
field due to small-scale nonlinear motions. Here we are not using 
linear theory to model deviations from a single global linear Hubble 
law, so the addition of velocity noise to our analysis is unnecessary. 
In our framework we would still have to take into account the noise 
due to peculiar velocities of galaxies within gravitationally bound 
clusters. However, in the COMPOSITE data set this has already 
been accounted for in gravitationally bound systems by assigning 
distances and associated uncertainties to clusters, rather than to 
the individual galaxies within the clusters. 

^ We use an overbar for uncertainties in the Hubble constant 
obtained by either radial or angular averages, to distinguish them 
from the distance uncertainties in individual data points. 



where = (Zlil'i ^) (Zlil'i ^) is the weighted 
mean distance of the Ng points in the 5th shell and 
cTq = 0.201 h~^Mpc is the distance uncertainty arising 
from the 20kms~^ uncertainty in the heliocentric pecu- 
liar veloc ity of both the Loc al Group and Local Sheet as 
given by iTullv et al.l (|2008f ) added in quadrature to the 
0.4% uncertainty in the mag nitude of the CMB dipol^E 
(|Fixsen et al.l ll996\ 

The uncertainty (j4]) is included since the Hubble law is 
necessarily determined by a linear fit through the origin 
for each shell. The local velocity uncertainty when di- 
vided by Hq provides an additional distance uncertainty 
in the mean distance of each shell relative to the origin, 
and the related uncertainty (|4]) in the mean slope Hg. 
This uncertainty is significant for shells close to the ori- 
gin, but much smaller for shells at large radii for which 
the mean distance has a long lever arm. 

In equations ^ and ^ Hq = 100 /ikms"^ Mpc"^ rep- 
resents the normalization used to convert velocity uncer- 
tainties to distance uncertainties in the COMPOSITE 
data set. One other important issue is the asymptotic 
value of the Hubble constant in each frame to which the 
variance in the Hubble fiow should be normalized, as this 
global value should have its own uncertainty. To this end 
we have divided the data in 12.5 h~^Mpc wide shells out 
to those radial distances of order 150/i~^Mpc, for two 
different choices of shells differing by the initial inner 
shell boundary, as shown in Table [Tl The penultimate 
shell, 10 or 10^ has been made wideo so that it contains 
a similar number of points to most inner shells. 

For both choices of shells 91 data points with r > 
156.25 /i~^Mpc in shell 11 have been used to determine 
the mean asymptotic value of the Hubble constant, H^, 
and its uncertainty. The inner boundary of this shell 
must be chosen at a sufficiently large distance that it is 
greater than the scale of statistical homogeneity. Thus 
we take the inner shell radius to be larger than the baryon 
acoustic oscillation (BAO) scale, that being the largest 
scale which could reasonably modify the gross features 

^ Very sl i ghtly diffe rent temperat u re di poles are given by 
Fixsen et aP ^99^ ) and Bennett et al.l ([200^ V Sin ce much of the 
COMPOSITE data set was determined before the IBennett et al.l 
(2002) result, we assume that it has been normalized relative to the 
heliocentric frame using the Fixsen et al. (1996) value of '^Qyij^ = 
371kms-i in a direction i = 264.14°, b = 48.26°, which is the 
standard used in the NASA/IPAC Extragalactic Database. Our 
heliocentric velocities of the LG and LS are taken from lTullv et all 
(|2008l ) as Vj^Q = 318.6 km s~^ towards £ = 106°, b = -6°, and 
^LS ~ 318.2 km s"-*^ towards i = 95°, b = —1° respectively. 

We originally included three shells in the range covered by shell 
10 or 10^ but have combined these at the request of the referee. 
This in no way changes any of our conclusions; if three shells are 
used then the uncertainties in Hs are simply somewhat larger on 
account of having fewer data points per shell. 
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of the local Hubble flo^. Furthermore, in the CMB 
frame the asymptotic Hubble constant should match 
the 100 h-^Mpc normalization used in the COMPOSITE 
data set. This is indeed satisfied by our choice. We find 
= (100.1 ± 1.7)/ikms-^Mpc"^ for the CMB frame 
and = (101.0 ± 1.7)/ikms-^ Mpc~^ for the LG and 
LS frames. We thus see that although the LG/LS value is 
1% larger than the CMB value, both values agree within 
uncertainties, and also with the distance normalization 
assumed in compiling the COMPOSITE sample. The 
fact that there are only 91 points in the outer shell is not 
an issue statistically, since the goodness of fit statistic is 
0.999, and angles on the sky are well sampled. 

We do not determine from the whole COMPOSITE 
sample, since it is dominated by points in the foreground, 
with a mean weighted distance of 15.05 h~^Mpc. The fit 
of a single linear Hubble law to the whole sample of 4,534 
points gives = (108.9 ± 1.5) /ikms~^ Mpc~^ in the 

CMB frame or = (104.4 ± 1.4) /ikms"^ Mpc"^ in the 
LG and LS frames. It is precisely because voids domi- 
nate the volume of space that we expect radial averages 
dominated by objects at sc ales less than t he diameters 
of the largest typical voids (jHovle fc Vogeley 2002, 2004; 
iPan e t al. 2011) to skew the simple linear average to val- 
ues greater than the asymptotic global value. This is 
confirmed by the full sample simple linear fit. Our pur- 
pose is to more carefully quantify the foreground Hubble 
fiow variance. 

The key statistical point about the determination of 
the mean asymptotic value, H^^ in each case is that its 
uncertainty provides a significant contribution to the to- 
tal uncertainty in the relative variation of the Hubble 
parameter in the 5th shell 



(5) 



We have checked that the angular sky coverage of the 
sample is consistent in individual shells. This is an im- 
portant check since we could get spurious results if the 
data in any shell was limited to one side of the sky, and 
potentially dominated by particular structures. In Fig- 
ure[T]we plot figures similar to Figure 1 of FWHIO, which 
shows the sky coverage within each of the unprimed shells 
of Section [2l We use a Mollweide projection in galactic 
coordinates (£, b) with £ = 360° on the extreme left and 
^ = 0° on the extreme right. Additional peculiar ve- 
locity information is encoded in the relative sizes and 
colours of the data points. We see that angular sky cov- 
erage is consistent in all almost shells, with some large 
gaps only the innermost shell 1, r < 12.5/i~^Mpc. For 
the primed shells, where the inner boundary is offset by 
6.25 h~^Mpc there is no sky coverage problem. We will 
retain the (unprimed) shell 1 in our analysis, but our 
statistical conclusions will not rely on it. 

^ Since the BAO feature is observed in galaxy clustering statistics 
in accord with the expectations of linear perturbation theory on 
a FLRW background, the scale of statistical homogeneity should 
necessarily be of the same order or smaller. In seeking convergence 
of bulk flows to the CMB dipole, researchers working in the peculiar 
velocity framework are currently considering the influence of the 
Shapley Supercluster on our local motion. Since Shapley is at a 
distance of 138 /i~^Mpc, such a large scale correlation would have 
to represent a very unusual fluctuation relative to the statistical 
homogeneity scale if the standard framework were correct. 



2.2. Results 

The results of our analysis are shown in Figure [21 where 
we plot SHs for both the CMB and LG frames, for the 
two independent choices of shells given in Table [TJ We 
computed the result for the LS also; it is essentially in- 
distinguishable from the LG frame. 

It is clear that the variance of the spherically averaged 
LG or LS frame Hubble flow is less than that of the 
CMB frame. In both frames the Hubble flow averaged in 
spherical shells gives SHg within 2a of uniform in almost 
all shell^ for fg > 55.1/i~^Mpc. However, particularly 
for values fg < 37.6/i~^Mpc, the LG/LS frame has SHs 
much closer to uniform than the CMB frame, and the 
average LG/LS frame flow is even within 1.2cr of uniform 
in the range 12.3 < ^5 < 23.5/i~^Mpc, whereas the 
average CMB frame flow is 4.0cr - 7.0a from uniform in 
the same range. Since the Local Sheet itself is deflned 
within r < 6.25 /i~^Mpc, this is not a result that would 
be readily expected with the standard interpretation of 
peculiar velocities. 

The statistical signiflcance of the relative uniformity of 
the averaged flow in the two frames can now be deter- 
mined by summing the mean square differences from a 
uniform 5H = expectation. 
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for each choice of rest frame in all shells outside an inner 
cutoff shell, Tg, as the inner cutoff is varied. An inner 
cutoff is commonly applied to eliminate the contribution 
of large peculiar velocities near the origin, and given the 
reinterpret ation we follow in this paper, the effect of vary- 
ing the cutoff is particularly interesting. The probability, 
^cmb(^s) ^lg(^«)' ^ uniform Hubble flow for each 
choice of rest frame and cutoff is then calculated directly 
from the complementary incomplete gamma function for 
the chi square distribution with the relevant number of 
degrees of freedom. A Bayes factor B{rs) = ^lg/^cmb 
for each choice of inner cutoff is determined for the two 
independent choices of shells in Table [TJ The resulting 
values of In B are tabulated in Table [TJ and plotted as a 
function of Vg in Figure [3l We also determined Pj^s/ -^cmb 
for the LS relative to the CMB; however, the values ob- 
tained gave Bayes factors which were essentially indis- 
tinguishable from those tabulated for the LG relative to 
the CMB. 

Figure [3] reveals a number of interesting features. 
The fact that overall the LG frame is more uniform 
is demonstrated by \nB being everywhere positive. If 
we consider a large inner cutoff, Ts > 106.25 /i-^Mpc 
then the difference in uniformity of the two frames has 
\nB < 1, which is not statistically signiflcant. With 
cutoffs in the range 37.5 < Vs < 100/i~^Mpc we 
flnd 1 < ln5 < 3 with positive evidence in favour of 
the LG frame being the more uniform. Bringing the 
cutoff down to Ts = 37.5/i~^Mpc gives \nB = 3.6 in- 

^ The one small exception is that ^Hq^^j^ is 2.1a from uni- 
form for shell 9 with 100 h'^ < r < 112.5 /i-^Mpc and fg = 
105.0 /i~-'^Mpc. In general, the LG frame flow is still somewhat 
closer to uniform than the CMB frame flow in the outer regions. 
For all shells with fs > 69.2 /i~-'^Mpc the LG frame flow is within 
1.36(7 of uniform. 
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creasing the Bayesian evidence to strong. For cutoffs 
Ts < 25 h~^Mpc the Bayesian evidence becomes very 
strong, liiB > 5. Different adjectives are used to de- 
scribe the stronge st Bayesian evidence (jKass fc Raftervl 
[19951 : iTrottal 120071 ): since \nB > 10 for any inner cutoff 
with Vs < 14.5 h~^Mpc Jeffrey's terminology of "deci- 
sive evidence" in favour of the relative uniformity of the 
Hubble flow in the LG frame seems to be appropriate. 

2.3. Nonlinearity from foreground structures and 
statistical issues 

We must be careful, however, in the determination of 
statistical confidence, since there are also statistically sig- 
nificant departures from uniformity in the LG and LS 
frames also, as is consistent with the presence of fore- 
ground structures. 

The nonhnear effects of the foreground structures can 
be seen by computing the goodness of fit probability, Qg, 
given by the complementary incomplete gamma function 
for Xs in shell s with u = Ng — 1 degrees of freedom. In 
Table [T] a bad linear fit is indicated in both the CMB and 
LG frames for shells s < 4 (unprimed) or 5 < 3' (primed) 
since Qg is less than 0.1 and equivalently the reduced 
per degree of freedom is significantly in excess of one. 

We have investigated the extent to which the relative 
magnitude of the Bayes factor is driven by the greater 
scatter relative to a linear law, rather than by the dif- 
ference of the linear fit of the Hubble constant from its 
asymptotic value. To investigate this in Tables [21 [3] and 
[4j we have repeated the analysis by successively removing 
points which contribute the greatest scatter: 

1. Firstly, in Table [2] we remove points which con- 
tribute an individual value Xi s = {'^i~^^i/Hs)'^ /erf 
with xfs > 5 in both the CMB and LG frames 
for both the primed and unprimed choices of shell 
boundaries. 

2. Secondly, in Table [3] we remove points with xf s > 5 
in either the CMB or LG frame for both choices of 
shell boundaries. 

3. Finally, in Table [4] we remove points with xf s > 5 
in either the CMB or LG frame for either choice of 
shell boundaries. 

Since the underlying Hubble relation is not linear, it 
is clear that we are rejecting some of the points with 
the strongest discriminating power in such a procedure. 
Indeed, only the points excluded in Table [2] could be said 
to be outliers in any sens43 since all other points have 
xf s < 5 in at least one frame and shell slicing. 

Nonetheless, although this procedure is not a perfect 
one it does illustrate that as the linear goodness of fit is 
improved the relative Bayes factor is somewhat reduced 
but still remains significant. In Table [2] only shell 4 has 
a significantly improved goodness of fit in both frames. 
However, in Table[3]the shells {2, 2', 3, 3'} now all have an 

^ The points excluded in Table [2] are still only outliers relative 
to the CMB and LG frames. Since we have not yet established 
a "minimum variance Hubble frame", it is perfectly possible that 
some of the points so excluded in fact have little scatter relative to 
a linear Hubble law in some other minimum variance frame. 



acceptable goodness of fit and a reduced x^ close to unity. 
Although In B is somewhat reduced, a very strong value 
ln5 > 5 is still found in shells {2, 2', 3} and indeed in 
shell 2 we still have \nB > 10. Our statistical conclusions 
are thus robust. 

In the final Table [4] even shells 1 and V have an accept- 
able goodness of fit. However, for shell 1 this comes at 
the expense of having removed two thirds of the 92 points 
originally present. The 12.5 h~^Mpc radius of shell 1 is 
simply so small relative to the foreground structures that 
a notion of a spherically averaged linear Hubble law is not 
really applicable. However, our statistical conclusions do 
not rely on the innermost shell. 

A better alternative for estimating Hs in the inner 
shells would be to use a robust optimization procedure. 
However, that would also require modeling the nonlinear- 
ity of the inner shells, and in the present paper we aim 
to be model-independent, without any particular theo- 
retical biases. To this end, we believe the very strong 
evidence for the relative uniformity of the LG and LS 
frames as compared to the CMB frame, has been clearly 
demonstrated. 

The outer shells with s > 5 (unprimed) or 5 > 4' 
(primed) all have a strong goodness of fit in the full data 
set of Table [TJ This is also true in the outermost, r > 
156.25 h~^Mpc shell, although it only contains 91 points. 
This outer shell, which is used to anchor the asymptotic 
value of the Hubble constant and its uncertainty, has an 
almost perfect goodness of fit Qs = 0.999 and a reduced 
X^ of 0.59 per degree of freedom in both frames. 

Some hints about the nature of the effects which con- 
tribute to the deviations from a uniform linear Hubble 
law can be discerned by comparing SHs in the shells 
where the deviations from uniformity become statisti- 
cally significant. Perhaps the most interesting feature is 
that counter to the overall trend, the individual CMB 
frame shells {4, 4', 5} with mean distances in the rang43 
41.0 < fs ^ 58.4/i~^Mpc have a closer to uniform 
Hs than the corresponding LG frame shells. In the cu- 
mulative Bayes factor this adds a negative contribution, 
and reduces the overall In 5 to 0.92 at fg =49.1 /i~^Mpc, 
even though adjacent points have \nB > 2. 

2.4. Systematic offsets from choice of reference frame 

Another important point is to consider how the non- 
linear dependence of Hs on the individual czi in the re- 
gression formula (j2j) can lead to systematic offsets when 
applying boosts. Suppose we are in a frame in which the 
spherically averaged variance in the Hubble flow is min- 
imized, which of course can be a frame other than the 
LG or LS one. Now change reference frame by applying 
a uniform boost to all data points, so that 

CZi CZ[ = CZi -\- V COS (7) 

where <pi is the angle on the sky between the data 
point and the boost direction. Then (cz^)^ i^^i)'^ = 
czf-\-2cZiV cos (l)i-\-v'^ cos^ in the numerator of (|2]), and 
cziri cz[ri = cziri + cos(^i) in the denominator. 

If we perform a spherically symmetric average ([2]) on 
data which is reasonably uniformly distributed over the 

Each bound is the average of the mean distances of the shell 
where the CMB frame is more uniform with the mean distance of 
the neighbouring shell where the LG frame is more uniform. 
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celestial spher^Hl then on average each positive contribu- 
tion of the linear terms v cos in the sums in the nu- 
merator and denominator of the boosted frame H'^ will 
be counterbalanced by a negative contribution from a 
V cos on the opposite side of the sky. The linear con- 
tributions are therefore roughly self-canceling, leaving a 
dominant contribution to the difference 



H[ — Ha 




{{v cos 



where {fi)^ 



(cZiTi) 



(8) 



)iY^i^i ^i ^) ^ is a weighted 



average. If we now consider (|8|) applied to successive 
shells with different values of 5, then given a uniformly 
symmetrical distribution of data the weighted average 
(^{v cos 1'^^ in the numerator will be roughly con- 

stant from shell to shell, while putting the leading order 
approximation czi H^Vj in the denominator we find 



(9) 



Consequently, for symmetrically distributed data the 
effect of incorporating a boost in the redshift data is 
an additional contribution to the spherically averaged 
Hubble constant which is inversely proportional to the 
averaged square distance. The difference between the 
CMB and LG frames in Fig. [T] does indeed show hints of 
such a dependence. Of course, the LG frame itself may 
incorporate such a dependence with respect to whatever 
frame has the minimum variance in i^^, only to a lesser 
extent. 

We stress that by our method of analysis the effect of 
a spurious boost is to add a spherically symmetric, or 
monopole, "Hubble bubble" type variation to the Hub- 
ble relation. This feature makes the present analysis 
very different to the standard peculiar velocity approach, 
where the focus is on dipole or higher multipole varia- 
tions. 

We summarize the results of this section as follows. Al- 
though there are significant foreground structures which 
distort the spherically averaged Hubble flow in a sta- 
tistically significant manner, the LG frame has a much 
smaller monopole Hubble flow variance than the CMB 
frame, counter to standard expectations. Nonethe- 
less, there is a particular range of distances at roughly 
40 ^ ^ ^ 60 /i~^Mpc for which the boost to the stan- 
dard CMB frame produces an apparently more uniform 
spherically symmetric average flow. This is the first ev- 
idence for the hypothesis we will present in Section [5l 
namely that rather than being a transformation which 
puts us in the frame in which the Hubble flow is most 
uniform at our own point, the boost to the CMB frame is 
actually compensating for the effect of foreground struc- 
tures largely associated with distance scales of order 

The absence of data in the zone of avoidance does not affect 
this argument, since the gaps in the data set are symmetrically 
distributed on opposite sides of the sky. The argument would only 
fail if there was a significant gap in the data on one side of the sky 
only. 



40 ^ to 60 ^Mpc. To better understand these struc- 
tures we now consider angular averages. 

3. ANGULAR AVERAGES 

In order to associate variance in the Hubble flow with 
particular foreground structures angular information is 
also required. The angular v ariance of the Hubble flow 
in the same HST K ey Data CFreedma n et al.ll2001[ ) in- 
vestigate d by iLi^ Schwarz (2008) has been studied by 
McCl ure fc Dved (|2007) (henceforth MD07) in the CMB 
reference frame. ^McCl ure fc Dyer (2007) used all 76 
points in the Key Data set, and concluded that a 13% 
variation in existed in the data. They based the sta- 
tistical significance of their conclusion on Monte Carlo 
simulations given the paucity of their data. Once again, 
although individual distances in the COMPOSITE sam- 
ple are noisier, the 60-fold increase in the size of data set 
enables a more detailed analysis. 

Since we wish to consider variations in by angle, 
we will firstly perform the most simple of radial separa- 
tions: we divide the data into an inner (r < Vq) and an 
outer (r > To) sphere, with a boundary Vq which we vary. 
We could simultaneously consider the angular variation 
within each of the radial shells of Section [21 given that 
the sky coverage is reasonably consistent from shell to 
shell, as seen in in Figure [H However, for a multipole 
analysis we will simply split the sample into two spheres, 
in order to determine the gross features of the relative 
angular variation with as much statistical confidence as 
possible. For the leading order dipole correction we will 
perform a split by radial shells in Section [X2l below. 

We will follow MD07 to produce contour maps of the 
angular variation in the average value of cz/r over the 
sky using a Gaussian window function. At each grid 
point on the sky, a mean is calculated in which the 
value of czi/vi for each data point is weighted according 
to its angular separation from the grid point. The Gaus- 
sian profile was given a standard deviation of = 25° 
by MD07, since with a narrower width the variation be- 
came dominated by the measurement errors, and angular 
resolution was compromised with a greater width. The 
weighting of the zth data point is given by 



1 



(10) 



where cos^^ = rgrid • = 25° is the smoothing scale 

and the Greek subscript a is used to represent the angu- 
lar dependence on the spherical polar coordinates, (^, h) 
encoded in 9i. 

We note from (p!Q|) that is the radius of a Gaussian 
window function. If the diameter 2(j^ is smaller than the 
angular size of the ZoA then windows centred on grid- 
points close to the galactic plane would have insufficient 
data to give reliable results in those regions. The ZoA is 
typically 30° wide for the COMPOSITE sample, which 
means that the smoothing scale must be greater than 
15°. On the other hand we cannot make the smoothing 
scale so large that we lose all angular resolution. This de- 
termines the choice = 25°, which matches that made 
by MD07. We have checked that varying the smoothing 
scale in the range 15° < < 40° does not significantly 
change the results that follow. 

Since (p!Q|) determines a mean value of Hq^ at each grid 
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point on the sky in which each data point is weighted 
by its distance from the grid point, there will be larger 
uncertainties for grid points near the ZoA, as can be seen 
in Figure [10] below. 

In the method adopted by MD07, the weighted mean 
Ho, is calculated at each spherical polar grid point by 



(11) 



with the weight (fTQ|) . The variance of this sample mean 
at each grid point is given by 



where 



CZiCFi 



(12) 



(13) 



is the standard uncertainty from error propagation of the 
uncertainty <Ji in the zth distance in (pT]) - (p!3|) pro- 
duced values. 

No additional uncertainty weighting of the weight Wi a. 
was used by MD07 in the determination of i^c^, since they 
claimed that the impact of the errors in the individual 
data points is averaged out by the Gaussian weighting 
procedure. The individual uncertainties in the COM- 
POSITE sample are larger, and therefore the question of 
the treatment of the uncertainty weightings in the deter- 
mination of the mean (pTj) is an important one. 

In order to manage the uncertainty weightings, rather 
than using equations (pT]) - (p!3|) . we will instead determine 
the weighted mean at each spherical polar grid point 
(i^ 6), by first evaluating its inverse 



(14) 



with the weight ([TO]) . The variance of ^ at each grid 
point is then given by 



where 



(15) 
(16) 



is the standard uncertainty in • = ri/{czi). Then 



(17) 



is the standard uncertainty in H^- If, following MD07, no 
additional uncertainty weightings are used then in ([14)) 
a is given by ()10p . Alternatively, if inverse variance 
(IV) uncertainty weightings are used then we replace ([TO]) 

by 

Wia= . exp(^). (18) 



H7 



The reason that it is preferable to work with is a 
consequence of the dominant uncertainties in the COM- 
POSITE sample being those associated with the distance 



measure, cr^. In the case of the radial shells we chose to 

minimize [cr^^{ri — cZi/H)] with respect to H for 
the same reason. The Gaussian window averaging adds 
a nonlinear weighting to what is otherwise a linear regres- 
sion. Using ([TS]) ensures that the nonlinear weighting is 
added to IV weightings determined from uncertainties 
([T6]) which are themselves linear in the measurement un- 
certainties. The alternative procedure of ([TT]) - ([T3l) intro- 
duces a different distance weighting of each point in ([T3]) 
which leads to different result^ when combined with the 
Gaussian window averaging in an equation analogous to 

We note that including a standard IV weighting in ([TS]) 
may not be the most robust method for uncertainty es- 
timates for this data set. In particular, as has been dis- 
cussed in WFH09 and FWHIO in the standard pecuHar 
velocity framework the nonideal geometry of typical sur- 
veys can lead to an aliasing of small scale power. Where 
the data is sparse biases can be introduced relative to the 
bulk flow of a regular volume that one is ideally interested 
in. To deal wi t h thes e issues Watkins et al. (2009) and 
iFeldman et al.l ()2010[ ) have developed a minimum vari- 
ance weighting method with respect to the leading pecu- 
liar velocity moment amplitudes (dipole, quadrupole and 
octupole). 

In the present paper we are moving beyond the pecu- 
liar velocity paradigm and so the methodology of WFII09 
and FWHIO cannot be directly applied. Nonetheless it is 
clear that an analogous angular minimum moment vari- 
ance weighting methodology should be developed to ro- 
bustly constrain angular variance of the Hubble flow. We 
leave such an approach to future work. In the present pa- 
per we will simply identify the principal angular features 
of the Hubble flow variance using the Gaussian window 
averaging method of MD07 by applying ([T4)) - ([T6)) either 
with no IV weighting via ([TO]) or alternatively with an 
IV weighting via ([T8]) . The difference between these two 
methods will turn out to be sufficiently small that it does 
n ot affect the main claim of Section 13. 4[ 

iMcClure fc Dverl (|2007[ ) performed Monte Carlo sim- 
ulations to assess the significance of the variation that 
they found. While a similar analysis has been under- 
taken by one of us (Smale 2012), in the present paper we 
will use the more direct method of plotting the contours 
of the uncertainty aa - as given by ([TTj) using either ([TO]) 
or ([T8|) - on the same map as the Hubble flow contours. 
Such a plot of the angular variation of the uncertainty 
contains much detailed information, uncertainties being 
larger in some angular regions rather than others. 

3.1. Angular variance 

In this subsection we will flrst obtain a qualitative un- 
derstanding of the angular variance of the Hubble flow by 
considering the variance of using ([TO]) with = 25° 
and no IV uncertainty weighting. In Figure [4] we plot the 
contour maps of angular Hubble flow variance produced 
using the whole data set in a single sphere, in both the 
CMB and LG rest frames. The map for the LS frame is 

We found that using equations ([TT]) - ([T3]) in place of ([l4)) - ([T6)) 
gives results which differ very little from each other if IV weightings 
are not used. However, once IV weightings are included using 
(|13|) gives values of Ha with a mean which is 10% lower than the 
mean values determined from the averages in spherical shells. 
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very similar to that of the LG frame. 

Next we produced a series of sky maps in the CMB, LG 
and LS frames with the data divided into two spheres, 
with the boundary between the two sphere varying in 
steps. We show a subset of these sky maps in Figures [5] 
and [6] for the CMB and LG rest frames. The LS frame 
sky maps are once again very similar to the LG frame 
ones. In each case we show the angular Hubble variance 
for data within an inner sphere, r < To, alongside the cor- 
responding map for the data in the outer r > Vq sphere, 
with To taking the values 12.5 20 4Qh~^ and 
mh-^ Mpc. 

The maps are of course not entirely independent, as 
there is overlap of data between the outer shells for small 
To and the inner shells of maps with larger Tq. The extent 
of overlap of sources, and their angular distribution, can 
be determined roughly from the numbers given in Table [T] 
and in Figure [H where points within individual shells are 
shown. We have deliberately chosen to work with maps 
which are not independent, so as to show how power in 
the dipole is transferred from the outer to inner sphere 
as To is varied. 

The first observation we make is that although both 
frames reveal a dipole structure, the nature of the dipole 
has important differences between the two frames. In 
the CMB frame the difference between the inner and 
outer spheres is not very strong. In the outer sphere 
the two poles migrate from being both in the northern 
hemisphere in the r > 12.5 /i~^Mpc map to both being 
close to the galactic equator in the r > 60 /i~^Mpc map, 
while the poles in the corresponding interior spheres be- 
come localized to the northern hemisphere. However, the 
strength of the dipole feature does not vary significantly 
between the inner and outer spheres, nor with the vari- 
ation of the boundary Vq between the inner and outer 
spheres. The fact that both poles are in the northern 
hemisphere in most of the CMB frame plots also means 
of course that the dipole is less strong than for example 
in the r > 12.5 /i~^Mpc LG frame map, for which the 
poles are closer to 180° apart. 

By contrast to the CMB frame, in the LG frame there 
is a significant radial dependence to the Hubble variance 
dipole evident in Figure [6l With the division set at Vq = 
12.5 /i~^Mpc there is very strong dipole feature in the 
outer r > 12.5 /i~^Mpc sphere, which is stronger than in 
the full sample map of Figure HI By contrast, within the 
inner 12.5 /i~^Mpc sphere any dipole signature is masked 
by other multipoles which appear equally as strong. 

As the division scale Vq is increased the relative power 
in the dipole in the inner sphere in the LG frame maps 
increases substantially, so that by the time we reach Vq = 
60 /i~^Mpc the inner sphere shows a dipole almost as dis- 
tinct as the outer sphere of the Vq = 12.5 /i~^Mpc map. 
At the same time the dipole in the outer r > 60 /i~^Mpc 
map becomes less distinct. This is consistent with our 
finding in the previous section that the structures princi- 
pally responsible for the Hubble flow variance lie within 
r < 65/i-^Mpc. 

Although there are less data in the outer shells, the 
angular sky coverage in the outer shells is very similar 
to that in the inner shells, as shown in Figure [TJ Thus 
the disappearance of a strong LG frame dipole in shells 
with r > 60 /i~^Mpc is not due to any spurious sampling 



issues. This conclusion will be confirmed by an indepen- 
dent analysis of the data in Section [3^ below. We remark 
that the dipole feature in the LG frame can be seen by eye 
in the colour coded peculiar velocities relative to Hs in 
each shell, as shown in Figure [T] in shells 2 and 3, which 
cover the range 12.bh~^ < r < 37.5/i~^Mpc there is a 
clear concentration of negative peculiar velocities (blue) 
in the upper left quadrant and positive peculiar velocities 
(red) in the lower right quadrant, which correlate with 
the dipole structure in Figure [6l These concentrations of 
peculiar velocities become more and more diluted by the 
contributions of peculiar velocities of the opposite sign 
in shells 4 and 5, where 37.5 < r < 62.5 h~^Mpc. In 
shells with r > 62.5 h~^Mpc the areas previously associ- 
ated with the dipole feature contain similar numbers of 
positive and negative peculiar velocities. 

The fact that the dipole in the CMB frame shows far 
less variation than that of the LG frame as Tq is varied is 
consistent with the hypothesis that it is not directly asso- 
ciated with the structures defining the Hubble flow vari- 
ance but is rather due to an overall systematic, namely 
the relative boost to the CMB frame, as discussed in 
Section [231 

The above statements are of course made from a sim- 
ple inspection of the sky maps by eye. However, the 
statements can be quantified by performing a spectral 
analysis on the sky maps. To this end we digitized 
the contour maps into 1 square degree regions and per- 
formed a multipole analysis using HEALPI}(0 to deter- 
mine the relative power in the Ci coefficients. On ac- 
count of the Gaussian window averaging there is aliasing 
at the 25° scale, and information for the high multipoles 
is not reliable. However, since multipoles with i > 4 
are very much suppressed a good measure of the sig- 
nificance of the dipole can be estimated by determining 
the quadrupole to dipole ratio, C2/C1, and octupole to 
dipole ratio, C3/C1, as listed in Table [5] in the inner and 
outer spheres as the boundary, Tq, is varied in the CMB, 
LG and LS frames. The inner and outer C2/C1 ratios are 
also illustrated graphically for the CMB and LG frames 
in Figure B 

In the LG frame C2/C1 = 0.061 in the outer r > 
12.5 h~^Mpc sphere, representing a small quadrupole rel- 
ative to dipole while C2/C1 = 0.653 in the corresponding 
inner sphere representing a quadrupole roughly compa- 
rable to the dipole. By contrast in the CMB frame with 
To = 12.5/i~^Mpc, the same ratio is C2/C1 = 0.096 in 
the outer sphere and C2/C1 = 0.123 in the inner sphere 
indicating a dipole which is similar in the inner and outer 
spheres, and less clearly defined than that of the outer 
LG frame. 

In the inner sphere the ratio C2/C1 in the LG frame 
drops substantially for Vq > 30/i~^Mpc, and maintains 
a value in the range 0.09 - 0.12 when AQh~^ ^ Tq < 
90 h~^Mpc. This is of course higher than the same ratio 
in the outer sphere; but the inner sphere value includes 
in every case a contribution from the innermost shell in 
which the dipole and quadrupole are comparable. 

While the dipole is generally less distinct in the CMB 
frame than in the LG frame, if we compare the specific 
values of C2/C1 in the other sphere when Vq = 40, 50 

http://healpix.jpl.nasa.gov/ 
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or 60 /i~^Mpc, we see that in each of these cases the 
dipole is stronger in the CMB frame than in the LG 
frame. The effect is strongest in the outer sphere for r > 
50 /i~^Mpc, when the quadrupole/dipole ratio drops to a 
minimum C2/C1 = 0.023 in the CMB frame as compared 
to C2/C1 = 0.032 in the LG frame. However, when To > 
70/i~^Mpc the dipole in the CMB outer shell becomes 
less distinct again. As we saw earlier the variance in the 
spherically averaged Hubble flow was less in the CMB 
frame in the range 35 ^ ^ ^ ^0 h~^Mpc. It appears 
that the boost to the CMB frame is also having the effect 
of making the angular variance of the CMB frame Hubble 
flow more dipole-like over this particular radial range. 

3.2. Constraining the dipole feature: Alternative 
methodology 

It is difficult to provide statistical bounds on the an- 
gular orientation and magnitude of the Hubble flow vari- 
ance dipole with Gaussian window averaging. However, 
a completely independent analysis can be made by fitting 
the raw data to a simple linear dipole law 



cz 



(19) 



for the LG and CMB rest frames, where in each case (f) is 
the angle on the sky between each galaxy and the direc- 
tion {id^bd) which defines the best fit dipole a xis. This 
method is similar to that use d in Figure 9 o f lSanda"g^ 
(fl986) or Figure 8 of Kashhnsi^it_aij (j2009b[ ). 

We note that more robust statistical results might be 
obtained if in place of (p!9|) we were to fit an alternative 
dipole law 

r 1 6^ 

- = — -^cos^', (20) 



CZ 



Hi 



since the uncertainties (p!6|) in — ri/{czi) are more 
directly related to measurement uncertainties than the 
uncertainties (p!3|) in Hi. However, while the two laws 
will agree if <C i^d, in general the relationship between 
(fT9|) and dSOj) is not linear, so that p ^ p' and ^ ^ 
when 13 /Hd can be of order 10%, as is typical for the 
data here. Since we have plotted the variation of Hq^ 
in Figures HHHl rather than ^ the only way we can 
expect to obtain angular agreement of the dipoles is to 
use (p!9|) . Furthermore, fitting (fT9]) should give results for 
the dipole in angular agreement with bulk flows found in 
the standard peculiar velocity framework. 

In each case we determine the four parameters Hd^ 
id and bd by a least squares fit of Hd and the linear pa- 
rameters Px = P COS id COS bd^ Py = P cos id sin bd^ and 
Pz = P sin bd. Details are given in the Appendix. We 
have performed the analysis in the CMB and LG frames 
in two ways: (i) in each of the radial shells defined in Ta- 
ble [TJ and (ii) with data aggregated into inner and outer 
shells split at a radius r^, as it is varied. In this way 
we have counterparts for both the analysis of section [2] 
and of section [3711 The results, by radial shell, are given 
in Table [6l and the magnitude of the associated dipole 
gradients is also plotted in Figure [H In Table [7] and Fig- 
ure [9] we show the equivalent results for all data outside 
a cutoff r > To as is varied. 

As in the case of the spherical averages the goodness 
of fit of the linear relation is poor in the first few radial 



shells, whose radius is smaller than the typical largest 
voids. We have checked that using the reduced datasets 
of Section [231 leads to a goodness of fit close to 1.0 with- 
out substantially changing any of the conclusions of the 
results that follow. 

We note that the values of Hd are smaller than those 
found in Section [2l This is a direct consequence of fitting 
the law (p!9|) , rather than the alternative (f2Q|) , and agrees 
with our observation in the case of the Gaussian window 
averages that when IV weightings are used a fit to (fTT]) - 
(p!3|) rather than (p!4|) - (p!6|) gives lower mean values of H^. 
We have checked that if (f2Q|) is used in place of (fig]) then 
the values of Hd agree within la with the values of Hg in 
Table [H for those shells with Qs > 0.2 in both fits. Thus 
it is the relative value P/Hd which is of most interest, 
rather than the absolute value of /3, in Tables [6] and [3 

First let us consider the results averaged outside r^. 
We see that in the LG frame there is a very strong 
dipole with magnitude p = 11.4 ± 0.4kms~^ Mpc~^ 
for To = 15/i~^Mpc, which decreases to a typical value 
P ^3.5 for the largest values of in Table [71 which of or- 
der 3<j different from zero. If we consider only those cases 
with goodness of fit Q > 0.1 then the most abrupt de- 
crease in P occurs in the range 40 < ro < 55 h~^Mpc 
where P decreases from 6.9 ± 0.4kms~^ Mpc~^ to 3.5 ± 
0.5kms~^ Mpc~^. The diminishing of the dipole feature 
at this scale in the LG frame is consistent with what 
was seen in the angular sky maps of the previous sec- 
tion. Furthermore the angular position of the dipole 
for 20 < To < 45 h~^Mpc is consistently in the range 
{id^bd) = (83° ±6°, -39° ±3°) while the dipole is strong, 
but the angular position then wanders once its magni- 
tude is reduced to residual levels. For > 80/i~^Mpc 
the typical position of the residual dipole differs from 
that of the inner dipole by 80° - 100° in galactic longi- 
tude. 

By contrast, the magnitude of P is initially smaller 
in the CMB frame for small values of To, with a value 
p = 6.1 ± 0.4kms-^Mpc"^ at 15/i"^Mpc which de- 
creases somewhat to P = 4.1 ± 0.5 at To = 40/i~^Mpc. 
However, P then increases to 5.6 =b 0.6 kms~^ Mpc~^ 
at To = 55 /i~^Mpc, and for larger its value re- 
mains roughly constant at this level, which is of order 
4<j - 7 a different from zero. Furthermore, the direc- 
tion {id^ bd) remains within la of the bulk flow direction 
{i,b) = (287°±9°,8°±6°) found by Watkins et al. (200i) 
for all values 20 < To < 115 h~^Mpc. For the largest 
values of in Table [7] the directiorO remains consistent 
wi th the bulk flow direct ion {i,b) = (319° ±18°, 7° ±14°) 
of iTurnbull et al.[ ([20121 ). 

In conclusion, a statistically significant persis- 
tent dipole feature remains in the outer spheres 
To >60/i~^Mpc in the CMB frame. In the standard 
peculiar velocity framework this dipole would be inter- 
preted as a bulk flow on the relevant scales. In the 
standard framework the whole sample is used, rather 
than splitting the sample into inner and outer spheres. 



14 In the outermost shell 11 of Table[6]with r > 156.25 /i-^Mpc, 
the uncertainty in the dipole position in the CMB frame is essen- 
tially the whole sky, meaning that there is not enough data to 
constrain the dipole in this range. In the LG frame the magnitude 
of the dipole in this shell is just 1.2a from (3 = 0. 
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The bulk flow is therefore usually determined as an 
average peculiar velocity magnitud e at the characteris- 
tic de p th of the samp l e; see , e.g.. iKocevs ki fc E belingI 
(j2QQ6[ ): lErdogdu et all (120061 ): iLavaux et a l. (2010). 

Our methodology is different here since we are not 
quoting a single bulk velocity magnitude at any particu- 
lar characteristic depth. For large Vq most of the data in 
the full COMPOSITE sample is not used, and the conse- 
quent uncertainties are reflected in the larger error bars 
for those distances in Figure [9l While there is less data 
in the outer regions, there is still enough data to give 
statistically significant results. 

Although a much less significant residual dipole re- 
mains for large Tq in the LG frame, the disappearance of 
the much larger dipole found at tq < 55/i~^Mpc in the 
LG frame indicates that this feature is associated with 
structures within this range. To pin down the distance 
scale of these structures, consider now the split by inde- 
pendent radial shells in Table [6] and Figure [HI In the two 
rest frames the magnitudes of the dipole agree in shell 3 
with = 30.2/i~^Mpc, and also in shell 5' with f^, = 
61.7 /i~^Mpc, but exhibit very different behaviour for the 
shells in between. In particular, the CMB dipole magni- 
tude reaches a minimum of /3 = 2.6 ± 0.6kms~^ Mpc~^ 
in shell 4, for which = 44.5 /i~^Mpc, whereas for the 

LG frame /3 = 14.9 ±0.8 km s~^ Mpc~^ in the same shell. 

We therefore have a statistically robust justifica- 
tion for the conclusion that the boost from the LG to 
CMB frame is compensating for structures in the range 
30 < r < 62 h~^Mpc. A large dipole in the LG frame 
is correlated with structures in this range, and the resid- 
ual dipole which remains in the LG frame beyond this 
scale is in a different direction. 

3.3. Dipole from Gaussian window averages 

Fitting a pure dipole (p!9|) to the Hubble flow leads to 
simple angular constraints on the dipole axis {£d^ ^d)- By 
contrast the representation of the uncertainties in the 
Gaussian window averaged sky maps is more compli- 
cated, since they represent a sum of several large angle 
multipoles. The uncertainties (p!7|) are most readily rep- 
resented as a function of angular position. In Figure [10] 
we show the example of the outer r > 15 h~^Mpc LG 
frame sky map, with the contours of the Hubble flow 
variance plotted as solid lines overlaid with colour map 
contours showing the angular uncertainties as a function 
of angular position on the sky. In each case we have used 
a, = 25°. 

The angular uncertainties are somewhat greater in a 
curved band near the galactic plane. This effect is due 
to the absence of data in the ZoA as well as the propa- 
gation of measurement uncertainties through the Gaus- 
sian window averaging procedure. We have checked that 
the same band of greater uncertainties is obtained in the 
CMB frame map, even though the positions of the ex- 
trema are quite different in that case. 

In Figure [To] we have identified the poles corresponding 
to the maximum and minimum of the Hubble variance 
in the LG frame, and have defined a la contour around 
each pole, where the la value has been taken as the max- 
imum on the map, i.e., 1.02/ikms~^ Mpc~^ in the un- 
weighted case and 0.75 /ikms~^ Mpc~^ in the weighted 
case. With the poles taken to be the positions of the 



red crosses in each case, we find (i^b) = (116°, —35°) 
and (i^b) = (249°, 21°) for the maximum and mini- 
mum respectively in the unweighted case. Similarly, 
{£, b) = (105°, -27°) and {£, b) = (253°, 24°) for the max- 
imum and minimum in the IV weighted case. We note 
that the poles are somewhat squeezed in galactic longi- 
tude as compared to a pure dipole. The latitude of the 
unweighted maximum matches that of the correspond- 
ing {id^ bd) = (68° ± 3, -38° ± 2°) from Table [7] within 
Icr, while the longitude of the minimum matches that of 
(£rf + 180°,-M. 

The mean value of obtained by Gaussian win- 
dow averaging, or alternatively by fitting a pure dipole 
as in Tables [6] and [T] differs in general from the val- 
ues obtained by spherical averages. For the r > 
15 h~^Mpc LG frame Gaussian window average, for ex- 
ample, is 97.3 /ikms~^ Mpc~^ in the unweighted case 

and 102.6 km s~^ Mpc~^ in the IV weighted case. The 
former value is close to the dipole averaged value of Ta- 
ble [T] and the latter to the spherical average. In the un- 
weighted case the maximum and minimum values of 
are +9.3% and -9.4% from the mean respectively. With 
IV weightings the differences are +11.0% and -7.9% re- 
spectively. In each case these differences are considerably 
larger than the standard uncertainty at any angle, which 
is of order 1% of the mean H^. The maximum varia- 
tion is comparable to the ratio P / Hd in Table [T] which 
is 12.0 ± 0.5% for r > 15 h-^Mpc. 

Finally, we remark that the strength of the residual 
dipole for r > 85/i~^Mpc might well change if higher 
multipoles are included. In particular, from Table [6] we 
see that in the outer independent shells the dipole mag- 
nitude in the LG frame is within 2(j of /3 = in shells 7^ 
8, 8^ 9, 10, 10' and 11. The average residual value in Fig- 
ure [9] in the LG frame is therefore largely driven by the 
data in shells 9^ 10. However, the goodness of fit in shell 
9' is only 0.40 for a dipole law as compared to 0.996 for 
the monopole fit of Table [TJ while Figure [3 show that the 
quadrupole to dipole ratio in the Gaussian window aver- 
ages begins to increase significantly at To = 100 h~^Mpc. 
The mean distance of shell 9' is 112/i~^Mpc. This is 
precisely the distance at which an effect might be ex- 
pected if the wall-to-wall distance-redshift is modified 
at the BAO scale, and if we are near the edge of a wall. 
A large increase in data at the relevant scales would be 
required to confirm this possibility. The magnitude of 
this effect is considerably smaller than that of the dipole 
at r < 55/i"^Mpc. 

3.4. Identification of Hubble variance with particular 
structures 

A large dipole structure in the Hubble flow across the 
sky is consistent with a foreground density gradient lead- 
ing to concentrations of more rapidly expanding void 
regions in one sector of the sky, and less rapidly ex- 
panding wall regions in the opposite sector. A detailed 
understanding of the structu res within 3 /i~^M pc may 
be gained from the work of iTullv et al.l (|2008l ) , and a 
viewing of the associated videc[3 is particularly instruc- 
tive. Sky maps of structures on larger scales are given by 

|http : / / if a. hawaii . edu / ^ tully / pecv-12min_sound-qt . mov| 
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lErdogdu e t J] (l2006l) using the 2 Micron All- Sky Red- 
shift Survey (2MASS). 

The distinctive feature of the location of our galaxy is 
that it is in a thin filamentary sheet, the LS which defines 
the supergalactic (X, F) plane, on the edge of a Local 
Void of at least 30/i~^Mpc diameter. While large void 
regions dominate one side of the sky, wall regions dom- 
inate the other side of the sky with the superclusters of 
Centaurus, Hydra and Norma being particularly promi- 
nent. Our Local Sheet and nearby filamentary sheets 
such as the Leo Spur are of modest density. The Virgo 
Cluster appears to be the closest region of the thick sec- 
tion of a dense nearby wall; however, it lies almost in the 
supergalactic plane of the LS, rather than along the axis 
which defines the greatest density contrast between the 
nearby voids and walls. 

In this particular situation our galaxy is neither in the 
middle of one of the largest typical voids of 30/i~^Mpc 
diameter, nor is it in the middle of one of the thick wall 
regions. Rather it in the transition zone between void 
and modest wall structures, close to the edge of both. 
In this circumstance the observed dipole pattern of Hub- 
ble flow variance might be expected to be the dominant 
one. Since the spatial width of typical walls is generally 
smaller than the diameter of the largest typical voids, ob- 
servers located in the middle of a thick wall region with 
extent in their (X, F) plane, with larger typical voids 
some way off and equidistant along their ±Z axes, might 
in contrast to our situation expect to see a more domi- 
nant quadrupole Hubble variance. 

The angular extent of various structures must also be 
important in determining how close the pattern of Hub- 
ble flow variance is to a dipole. A pure step function 
contains many higher multipoles, so a simple division of 
the whole sky into two hemispheres of uniform faster and 
slower expansion would contain many higher multipoles. 
Since voids have a purer ellipsoidal geometry than walls, 
in terms of defining the relevant angular scales it is the 
voids which will more clearly delineate the dipole density 
gradient. 

Some estimate of the angular scales of the neare st voids 
can be obtained from the work of iTully et al.l (j2008). 
However, since the dis tances of ni a ny ga laxies in their 
survey are not known iTullv et all ()2QQ8l ) present their 
diagrams in redshift space, which are subject to redshift 
space distortions as large as the Hubble flow variance 
that we are endeavouring to measure. 

Our own Local Void comprises three separate smaller 
volumes: the Inner Local Void, and the Local Voids 
North and South, w hich are separated by filamentary 
thin wah structures (jTullv et al.l I2QQ8D . Here "north" 
and "south" refer to directions orthogonal to the LG 
in supergalactic coordinates; and since the plane of our 
galaxy is roughly orthogonal to that of the LS, this means 
that supergalactic "north" and "south" indicate direc- 
tion principally along the galactic longitude axis relative 
to the supergalactic north pole at £ = 47.37°, b = 6.32°. 

The Inner Local Void, which is ellipsoidal with its ma- 
jor axis roughly parallel to the Local Sheet, is the struc- 
ture that covers the largest fraction of the sky in the Lo- 
cal Void complex. From Figure 10 of (jTullv et al.ll2QQ8[ ) 
we estimate that it covers at least 40-60% of one hemi- 
sphere, given the uncertainties of redshift space distor- 



tions. In any case we expect it to be too large a fraction 
of the sky to give a pure dipole. This is confirmed by 
splitting the inner and outer shells at Tq = 12.5 /i~^Mpc, 
since the inner shell should just exclude the Inner Local 
Void while retaining the Local Voids North and South. 
As shown in the first panel of Figure [6] and in Table O 
the inner sphere in these cases has similar power in the 
quadrupole. 

The dipole axis appears to be principally defined 
by structures within the range 30 /i~^-62 /i~^Mpc, 
which lies beyond the scales explicitly identified by 
ITullv et all (I260I However, using the 2MASS sur- 
vev lErdogdu et all (|2006[ ) have reconstructed the density 
field in shells every 20/i"^Mpc out to 160/i"^Mpc. To 
define a dipole, rather than simply locating the largest 
overdensity or underdensity, one must find an axis where 
the integrated density gradient, including foregrounds, is 
maximized. If we compare the results of Sections 13.31 to 
Figure 3 of lErdogdu et al.l (2006) we see that the mini- 
mum Hubble variance pole coincides with the near side 
of the Centaurus-Hydra Wall on one side of the sky and 
the maximum Hubble variance pole coincides with the 
Andromeda Void on the opposite side of the sky. Our 
axis to the Andromeda Void passes through the Inner 
Local Void and the edge of the Local Void NortlT^. 

The centre of the Hydra Supercluster is at r = (34.9 ± 
2.5)/i-^Mpc, {£,b) = (269.6°, 26.5°), and the centre of 
the Centaurus Supercluster is at (31.5 ± 2.6) /i~^Mpc, 
{£,b) = (302.4°, 21.6°). These radial distances are close 
to the radial scale at which the CMB frame Hubble 
flow begins to be more unifor m than the LG one . In 
the next sky map plotted by Erdogdu et alj (|2006l ). at 
r = 40/i~^Mpc, the Hydra Supercluster remains very 
dense near the Hubble variance angular minimum while 
on the opposite side of the sky the Andromeda void 
has begun to close up, but an adjacent void complex 
Cygnus- Aquarius has opened up, maintaining the dipole 
density gradient. By Figure 5 of Erdogdu et al. (20Q6|) 
at r = 60/i~^Mpc, on the other hand, there are now 
large over densities, Pegasus and Pisces, in the angu- 
lar patch that previously contained the Hubble variance 
maximum, while on the opposite side of the sky in Hydra 
underdensities have emerged. These opposing influences 
will now even out the density gradients along the dipole 
axis. We can therefore understand why the dipole dimin- 
ishes beyond r = 60 h~^Mpc. 

4. CORRELATION OF HUBBLE VARIANCE AND CMB 
DIPOLES 

Having demonstrated that Hubble flow is more uni- 
form in the LG and LS frames as compared to the CMB 
frame, and that there is a persistent dipole feature in the 
variance of the Hubble flow when viewed in the LG and 
LS frames, the natural question to ask is to what degree 
is the Hubble flow variance dipole correlated with the 
CMB dipole? 

16 In the terminology of ITullv et aP f200F) the "Local Void 
North" comprises the reg i on de noted the "Delphinus Void" in Fig- 
ure 3 of Erdogd u et all (|2006l ) together with an adjacent large 
6 < area extending to just aloove the galactic plane, 6~6, with 
47 < £ < 90. The "Local Void South" of T ullv et al.l (|2008i ) is 
similarly much larger than the area marked "LV" in Figure 3 of 
Erdogdu et al. (2006) and extends to adjacent 6 < regions above 
the galactic plane, with i < 47. 



12 



Wiltshire, Smale, Mattsson & Watkins 



To answer this question we must compensate for our 
hehocentric motion with respect to the rest frame of the 
LG or LS by performing a boost to the relevant rest frame 
and examine the residual CMB temperature dipole in the 
rest frame in question. We create an artificial residual 
CMB dipole temperature map by subtracting a boosted 
CMB sky with temperature 



T 



7(1 - (v/c) COS0') 



(21) 



from the corresponding observed pure temperature 
monopole plus dipole maps using the iFixsen et alJ 
(|T996I) values assumed in Section [2l Here v = v^^ 
OT V = v^^ as appropriate, and 7 ~ 1 since ve- 
locities are nonrelativistic. This leaves us with a 
residual temperature dipole with poles ±5.77 mK at 
{i,b) = {(97.0, -29.3), (277.0, 29.3)} in the LG frame, 
and ±5.73 mK at {i, h) = {(90.0, -26.6), (270.0, 26.6)} in 
the LS frame. The dipole amplitudes have a 6.3% un- 
certainty arising principally from the uncertainty in the 
heliocentric velocity of the LG and LS frames. The LG 
residual temperature dipole is shown in Figure [TTJ 

We compute a correlation function directly by us- 
ing HEALPIX to digitize both the residual temperature 
dipole map, and also the corresponding Hubble fiow vari- 
ance maps for the LG or LS frame as relevant. We quan- 
tify the correlation between the variance of that Hubble 
expansion and the residual CMB temperature dipole by 
the Pearson correlation coefficient 
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(22) 

where is the temperature in the pixel with angular 
coordinates a, T is the mean temperature. 
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Ha is given by ([11]), aa by (fT7|) . and Np denotes the total 
number of pixels distributed over the sky. As we are con- 
sidering a pure residual CMB temperature dipole there 
are no uncertainties in T^. Since HEALPIX partitions 
the celestial sphere into pixels of equal area, and since 
the CMB temperature dipole is assumed to be ideal, the 
only weighting in the sum comes from the measurement 
uncertainties of the Hubble ffow. 

With (jg = 25° we performed a correlation analysis be- 
tween the Hubble variance dipole and the residual CMB 
temperature dipole in both the LG and LS frames, as the 
division radius, Tq, between the inner and outer spheres 
was varied. The results are shown in Figure [121 We 
observe firstly that the correlation coefficient is negative 
since the maximum value of the Hubble parameter co- 
incides with the minimum residual CMB temperature. 
The strongest anticorrelation is therefore represented by 
those values which are closest to —1. 

In all cases the result for the LS frame does not differ 
greatly from that of the LG frame. However, the anticor- 
relation is generally a bit stronger in the LG frame. The 



anticorrelation is stronger for the IV weighted sky maps 
in both frames. The anticorrelation is strongest in the 
outer sphere for To = 15 h~^Mpc. As might be expected 
from Section 13. If the anticorrelation remains strong in 
the outer sky maps for values of To up to 40 h~^Mpc. By 
contrast, the anticorrelation in the inner sphere is not 
at all strong for r < Tq with small values of Tq. How- 
ever, the anticorrelation in the inner sphere improves 
dramatically as Tq is increased, and by the stage that we 
reach Tq = 50/i~^Mpc the anticorrelation is comparable 
in both spheres for the unweighted case, and stronger 
in the inner sphere than in the outer sphere for the IV 
weighted case. There is no further improvement in the 
anticorrelation in the inner sphere for Tq > 60/i~^Mpc, 
which is again consistent with the earlier indications that 
the structures responsible for the Hubble variance dipole 
are within 65 h~^Mpc. 

One final question is the extent to which the correlation 
depends on the Gaussian smoothing width, a^. We have 
checked this is two ways. Firstly, we have recomputed the 
correlation coefficient for a range of values of for the 
r > 15 h~^Mpc map, the case which shows the strongest 
anticorrelation. The results are shown in Table [51 We 
find that the anticorrelation in the IV weighted map is 
stronger than the unweighted map for different choices 
of (7^. Moreover, as well as the correlation coefficient 
in the LS frame being slightly weaker, it also varies a 
little more with smoothing angle. In the LG frame the 
correlation coefficient varies the least for varying in 
the IV weighted case. Indeed to two significant figures 
the correlation coefficient of —0.92 is unchanged for = 
25° ±5°. 

A second check on the relation between the CMB tem- 
perature dipole and Hubble ffow variance, that is com- 
pletely independent of the Gaussian window averaging 
procedure, is given by evaluating a correlation coefficient 
between the residual temperature dipole and the raw 
COMPOSITE peculiar velocity data in the LG frame. 
In this case the relevant correlation coefficient is given 

by 
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where Vi denote the peculiar velocities and TV is the num- 
ber of data points. The weighted average peculiar ve- 
locity should approach zero for a large number of data 
points: here v = — 64.9kms~^ with a standard deviation 
of 722.4 km s-^ For the COMPOSITE LG-frame veloci- 
ties from the data with r > 15 /i~^Mpc, N = 4359 and we 
obtain = —0.35. The magnitude of this correlation 
is naturally lower than it is for the weighted grid data 
which we calculated above due to the scatter in these 
data, but the well defined number of points implies we 
have better statistical tools to quantify our confidence 
that the correlation is indeed nonzero. We test this by 
evaluating the variable 
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where v = N — 2 is the number of degrees of freedom. If 
there is no correlation, the test variable t in (j25j) should 
follow the standard normal distribution A^(0, 1). For u = 
4357 and = —0.346, we obtain the value t = —24.35, 
i.e., a deviation of more than 24cr. This is extremely 
strong statistical evidence for a nonzero (anti) correlation. 

We therefore have strong evidence that the dipole fea- 
tures of the two maps in Figure [TT] are related. From Fig- 
ure [10] we see that in the IV weighted sky map the cooler 
residual CMB temperature pole (marked by a cross) lies 
just 7.4° from the maximum of the Hubble flow variance, 
well within the la contour. The hotter residual CMB 
temperature pole is separated by 22° from the minimum 
of the Hubble flow variance, however, and lies 10° outside 
the 1(7 contour but is within 3cr. It remains to be seen 
whether the uncertainty estimates in this case are signifi- 
cantly affected by the choice of weighting scheme. In par- 
ticular, the left hand panel of Figure 1 of FWHIO shows 
that with a maximum likelihood estimate based on tra- 
ditional IV weightings there are a substantial number of 
very strongly weighted data points in the COMPOSITE 
sample in the region which coincides with that of mini- 
mum aa to the north of residual CMB temperature pole 
in Figure [TOT b). With the alternative minimum variance 
weightings shown in the right hand panel of Figure 1 of 
FWHIO the same data points are not strongly weighted. 

5. ORIGIN OF THE CMB DIPOLE 

Ever since the first bounds were placed on 
the anisotropy of the CM B in the mid 1960s 
dPartridge & Wilkinsojil 11967^ it has been assumed 
that the dipole anisotropy represents a measurement of 
our motion with respect to our surfa ce of average ho- 
mogeneity (jPeebles fc WilkinsonI 119681 ). Yet our results 
clearly show that for spherical (monopole) averages the 
Hubble flow is closer to being uniform in the frame 
of the LG or LS, rather than the frame indicated by 
treating the CMB dipole as being entirely due to a 
boost. This is completely unexpected in the standard 
framework, since the cosmic rest frame and the frame of 
minimum Hubble flow variance should be one and the 
same. Moreover, we find a strong dipole in the variance 
of the Hubble flow correlated with the residual CMB 
temperature dipole. 

If the standard peculiar velocity framework were cor- 
rect, then conceivably the residual bulk flows might give 
a dipole pattern in the Hubble flow variance consistent 
with our observation of the dipole variance. Indeed, in 
the peculiar velocity frar nework th e dipo le pattern we 
observe was detected by iJha et all (|2007[ ) in a sample 
of 69 MLCS2k2-reduced Snela distances in the range 
0.005 < z < 0.025, and interpreted by them as a pos- 
itive detection of the motion of the LG with magnitude 
541±75kms-i towards {£,b) = (258° ± 18°, 51° ± 12°), 
which is 2<j consistent with the amplitude and magni- 
tude of the boost of the LG with respec t to the CMB 
frame. The problem was re-examined by iGordon et al.l 
using a sample of 61 SALT-reduced Snela in the 
range 0.0076 < z < 0.124, who found v^^ = 697 ± 137 
in the direction {£,b) = (220° ± 14°, 32° ± 11°). When 
the linear theory was corrected to account for corre- 
lated peculiar velocities, the uncertainties on these val- 
ues were increased giving v^^ = 690 ±201 towards 
{i,b) = (257° ± 24°, 29° ± 16°) ([Gordon et al.l [20081 



which is now la consistent with the LG boost with re- 
spect to the CMB frame. 

Phenomenologically, the above results in the linear the- 
ory sho w agreement with the expected LG boost, and the 
iJha et al. (2007) Snela sample appears to show conver- 
gence to the expected result within the same radial scale 
we find in the COMPOSITE sample: in Section [2] we saw 
that the LG frame is within 1.36<j of uniform in spherical 
shells with fg > 69/i~^Mpc, and in Section [321 we saw 
that the large dipole feature in the LG frame diminished 
to its residual value by To = 55 h~^Mpc. Thus the struc- 
tures responsible for both the monopole and the dipole 
variation in the Hubble flow appear to be foregrounds 
within roughl£![ 65 /^-^Mpc, or z 0.022. 

In the standard peculiar velocity framework, however, 
demonstrating the convergence of bulk flows on this scale 
has proved challenging. In particular, in the linearly per- 
turbed FLRW model the peculiar velocities are given by 
(lPeebleslll993[) 
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where (5^(r) = {p — p) / p is the density contrast, and the 
power 0.55 of the matter density parameter ^^q, rather 
than 0. 6, ^ives impro ved accuracy for models with dark 
energy (|Linder|[2Q05[ ). 

The search for convergence of bulk flows within this 
framew ork has a thre e deca de histo ry, which is summa - 
rized bv lLavaux etaU (j2010[ ) and b ylBilicki etaU (I2Q1JJ) - 
Contrary to earlier investigations (|Erdogdu et al. |2006[ ) 
iLavaux et al] (|2010 ^ failed to find convergence in the 
2MASS survey on scales up to 120/i~^Mpc: less than 
half the amplitude was generated on scales 40 /i~^Mpc, 
and whereas most of the amplitude was gen erated within 
120 /z~ ^Mpc the direction did not agree. iBilicki et al.l 
(|2011[ ) analyse a larger sample in the 2MASS survey us- 
ing a different methodology and fail to find convergence 
within 150 /i"^Mpc. 

In the COMPOSITE sample iWatkins et al[ (|2QQ9[ ) 
failed to see convergence of peculiar velocities to the 
CMB dipole on scales of 50 h~^Mpc. Even more puzzling 
in their results is the suggestion that the bulk flow actu- 
ally increases with increasing scale above 20 h~^Mpc (see 
their Figure 5). In the peculiar velocity framework, the 
only way to understand how a larger volume can have a 
larger bulk flow than a smaller volume contained within 
it is to posit that the inner volume has an additional 
compensating motion in the opposite direction. While 
possible, this arrangement seems unexpected at best. 

Given these problems and puzzles of the standard lin- 
ear theory, combined with our quite unexpected results 
from Sections [SHU the natural course of action is to 
go back to first principles. From the point of general 
relativity the framework that is commonly adopted in 
determining peculiar velocities has no obvious justifica- 
tion. In particular, it is well known that velocities are 
only strictly defined in terms of local Lorentz boosts at 
a point, and the extension of the concept to larger scales 

There is some ambiguity in defining this transition scale. Con- 
servatively, we take the scale to be the average of the mean dis- 
tances of shells 5' and 6, where 6Hs first drops below 1.5cr difference 
from uniform. 
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requires geometrical assumptions. In particular, the as- 
sumption that a local inertial frame can be extended to 
scales of tens of megaparsecs is arbitrary. Every known 
inhomogeneous exact solution of the Einstein equations 
exhibits differential expansion of space. A differential 
expansion of space should therefore be expected to be 
the natural state of affairs on scales which are deemed to 
contain nonlinear inhomogeneities in the standard cos- 
mological framework. 

Although there are many sources of CMB anisotropics, 
which are extremely well studied, there is one source of 
CMB anisotropy in a universe which is not considered - 
at least not in the form we propose. If the universe is 
assumed to have been homogeneous at the epoch of last 
scattering, with a close to uniform temperature at the 
time of photon decoupling, then temperature differences 
will be generated if the universe has grown sufficiently 
inhomogeneous by the present epoch that the distance 
to the surface of photon decoupling is slightly different 
in different directions. Over large distances photon paths 
see an average of all the structures, but the last section 
of a photon journey below the scale of statistical homo- 
geneity will be influenced by the peculiar foregrounds. 

Our universe is in fact sufficiently inhomogeneous on 
scales < 100/i~^Mpc that the differential expansion of 
void and wall regions can be expected to produce dif- 
ferences in the distance to the surface of last scattering. 
The important question of how large such an effect could 
be expected to be must necessarily involve cosmological 
model assumptions. Here we will not attempt a model 
dependent theoretical estimate; but rather will determine 
the magnitude of the effect required to produce the ob- 
served residual temperature dipole. 

As long as the average evolution of the universe can 
be described by an average cosmic scale factoiH which is 
related to the observed cosmological redshift by a^j a = 
1 + and which is inversely proportional to the mean 
CMB temperature, T oc 1/a, then a small change, Sz, in 
the redshift of the surface of photon decoupling - due to 
foreground structures - will induce a CMB temperature 
increment T = T^+^T, where = 2.725 K is the present 
epoch mean temperature and 
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Zdec = 1089 being the mean redshift of photon de- 
coupHng. For the LG residual dipole the increment 
6T = ±(5.77 ± 0.36) mK represents a redshift increment 
(5z = t(2.31±0.15). 

The comoving distance of the surface of photon decou- 
pling is given by 
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in the standard spatially flat ACDM model, where Qj^^ = 
1 - - ^i.n aud 1^ r>n = 4.15/i~^ X 10~^. If wc takc 
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It is not necessary for the average evolution to obey the Fried- 
mann equation for this statement to be true. In particular, it also 
applies to the timescape cosmology ( Wiltshire 2007a, b, 2009) which 
describes the average evolution of an ensemble of void and wall re- 
gions in the Buchert averaging scheme (Buchert 2000, 2008), with 
additional physical assumptions about the interpretation of physi- 
cal observables. 



^Mo ~ ^'^^ ^ ~ ^•'^^ ^^^^ comoving dis- 

tance increment of SD = ^(0.33 ± 0.02) h~^Mpc is what 
is required to generate the residual CMB dipole in the LG 
frame. For ^j^q = 0.30 the value is slightly reduced to 
SD = =f(0.32 ± 0.02) /i-^Mpc. For the timescape model 
(| Wiltshire! l2QQ9f ) the value is similar, with possible small 
differences depending on parameter values. 

The results of the previous sections suggest that the 
structures responsible for the Hubble flow variance dipole 
lie at most within 65 h~^Mpc. Since the differences in the 
distance to the surface of last scattering occur effectively 
at z = 0, a 0.35 h~^Mpc difference in distance therefore 
would amount to a maximum 0.5% difference on these 
scales. Even if the whole difference was taken up within 
an average distance scale of 30 /i~^Mpc, leading to a 1% 
effect, this is still within the regime of plausibility given 
the degree of Hubble flow variance we observe in the 
COMPOSITE sample. 

Our picture, therefore, is one of differential expansion 
of void and wall regions at late epochs leading to dis- 
tance differences of up to the order of 1% between walls 
and voids on 30 - 62/i~^Mpc average distance scales. 
While such differences are not isolated to our own imme- 
diate vicinity, when light travels over scales larger than 
100 h~^Mpc the differences generally average out on any 
typical line of sight. It is on the last stretch of the journey 
on average distance scales < 65 h~^Mpc that the partic- 
ular foreground inhomogeneities peculiar to our own lo- 
cation give a net anisotropic contribution. In our case the 
largest foreground density gradient defines an axis with 
the void direction yielding a slightly larger distance than 
average and a net CMB temperature decrement, and the 
opposite wall direction a slightly smaller distance than 
average and a CMB temperature increment. 

We should stress that an alternative origin for part of 
the CMB dipole which we propose here will in no way 
affect the interpretation of observations of the angular 
power spectrum on scales smaller than 15°, as the cor- 
relations that we are considering only involve larger an- 
gles. In particular, although measurements of the acous- 
tic pe aks have now reached high precision (jKeisler et al.l 
'2011) such observations deal with scales of less than 1°. 

One important question remains: why is the CMB 
dipole so large compared to other multipoles if a substan- 
tial contribution is due to an anisotropy in the distance- 
redshift relation arising from foreground structures? We 
will not attempt to answer this question here. The prob- 
lem of describing the propagation of light through a real- 
istic inhomogeneous structure is a complicated one, and 
in the case of CMB photons we are dealing not with sin- 
gle photons but a statistical distribution of them. All 
that we can say is that we have strong evidence for a 
monopole and dipole variance of the Hubble flow, and 
that the dipole variance is closely correlated with the 
CMB temperature dipole. 

Past investigations do support the plausibility of our 
proposal, however. Studies of ray tracing of the CMB 
sky as seen by an off-centre observer in a large Lemaitre- 
Tolman-Bondi void have been shown give to very small 
quadrupole and octupole coefficients a^^ and a^^ as com- 
pared to the the dipole coefficient a.^. For example, a 
Newtonian approximation numerically close to the ray- 
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traced result yields (|Alnes fc Amarzguiouill2QQ6l ) 

So _ / 15 {hin — ^out)<^off /^gx 

~ V 4 2998 Mpc ^ ^ 

where i^ino = 100 /iin kms~^ Mpc~^ and i^outo = 
100 /lout kms~^ Mpc~^ are the Hubble constants in- 
side and outside the void, and doff is the distance of 
the observer from the centre in Mpc. Although we 
are not d ealing with large voids here, the results of 
lAlnes fc Amarz guioui (2006) can still be applied. Even 
for the relatively large values doff = 50/i~^Mpc and 
hin — hont = 0.2, we find that ratio of <^2o/^io ^^^^ 
of order 1%. 

The actual matter distribution is of course much more 
complicated than that of LTB voids, and the problem 
of determining the average propagation of CMB photons 
through the foregrounds is closely related to how to real- 
istically average foreground density fields. Nonetheless, 
the simple LTB example shows foreground structures can 
generate a CMB dipole without simultaneously generat- 
ing other multipoles larger than those seen observation- 
ally. 

6. DISCUSSION 

In summary, we have shown with decisive Bayesian 
evidence that when averaged in spherical shells on scales 
< 150/i~^Mpc the Hubble fiow is more uniform in the 
rest frame of the LG or LS than in the standard "rest 
frame" of the CMB. This result is independent of any 
cosmological model assumptions other than the existence 
of a suitably averaged linear Hubble law. Next we found 
that the variance of the Hubble fiow has a distinct dipole 
structure in the LG and LS frames, especially when the 
closest foregrounds on scales < 15 h~^Mpc are omitted. 

Various pieces of evidence, including in particular 
the analysis of Section 13. 2[ suggest that the dipole 
originates from a peculiar density gradient on scales 
30 < r < 62 h~^Mpc. Since light traveling to us from 
distant regions of the universe must necessarily travel 
through the foreground structures, these structures will 
also lead to systematic corrections to the distance- 
redshift relation even on the largest scales. The corre- 
lation between the Hubble variance dipole and the resid- 
ual CMB temperature dipole suggests that the principal 
component of the CMB dipole, apart from a local boost, 
is indeed such a systematic. If the CMB dipole is inter- 
preted in this fashion then we can expect that on large 
scales ^ 100/i~^Mpc the typical distance for a given 
redshift might vary from the average distance by up to 
0.35 h~^Mpc. This result is likely to have important con- 
sequences for cosmology. In this section we will discuss 
related issues, which are necessarily speculative, but sug- 
gest many avenues for future research. 

The most profound conceptual issue is that we should 
be careful to think about the differential expansion of 
space that results from the expansion history of regions 
of different density. By tradition, redshifts on scales up 
to z = 0.1 are treated by a special relativistic Doppler law 
in Euclidean space, and variations in the Hubble fiow are 
seen as peculiar velocities with respect to this expansion. 
Redshift space distortions resulting from the variance 
of the Hubble fiow in the peculiar velocity framework 



are w ell understood in terms of the Kaiser effect (jKaiserl 
Il987l ). However, the peculiar velocit ies derived from (j26j) 
in the standard linear framework (jPeeblesI Il993f ) have 
not been found to converge to the observ ed bulk fiows 
within distance scales surveyed to date (|Lavaux et al.l 
2010; Bilicki et al.|l201lD. 

From general relativity there is no a priori reason for 
assuming that space is Euclidean with a simple Doppler 
law on scales < 150/i~^Mpc. Nonetheless this assump- 
tion is so firmly embedded in much of the practice of 
observational cosmology that it is nontrivial to disentan- 
gle the consequences of revisiting this assumption. Here 
we will sketch just some of the directions which should 
be pursued in more detailed investigations. 

6.1. From bulk flows to Hubble flow variance 

From a formal point of view, in order to fully imple- 
ment the Hubble fiow variance paradigm a new multipole 
expansion needs to be developed to statistically constrain 
the Hubble fiow variance. We will leave this for future 
work, but make a few preliminary observations. 

A dipole anisotropy of order ±0.5% in the distance- 
redshift relation on 65 h~^Mpc scales is of course con- 
siderably less than the typical values of P/ ^^9% for 
individual shells in the range 30 h~^ ^ 62/i~^Mpc, 
as given in Table [6l It is also considerably less than 
the 4-8% variation of P/H^ seen in the averages with 
r > To in the range 30 h~^ < Tq < 50/i~^Mpc as given 
in Table [71 or similar angular variations of the Hubble 
fiow discussed in Section [331 However, the values of the 
dipole magnitude in Table [71 are obtained by applying a 
simple dipole law to the whole sample, and similarly in 
Section 13.11 we considered the angular deviation from a 
single linear Hubble law on the whole sample. As we see 
in Figure [2fb) there is considerable radial variation in the 
averaged Hubble law, and ideally this must be factored 
out before additional angular variation is considered. 

In particular, the 0.5% dipole anisotropy in the 
distance-redshift relation on 65 h~^M.Y)C scales repre- 
sents a departure from the spherically averaged distance- 
redshift relation. The starting point for a multipole ex- 
pansion of the Hubble fiow variance is to define the aver- 
age comoving distance, I), to a redshift, z < ^hom, within 
the scale of statistical homogeneity, ^hom, according to 




Here Hs{zs) is computed exactly as in (|2]) except that 
the shells on which the linear regression is performed 
are chosen by redshift ranges, Zs < z < Zs ^ cf^^ 
where cr^ represents the width of the radial shells in 
redshift. For example, in Section [2l we chose shells of 
radial width 12.5/i~^Mpc, which corresponds to taking 
c(Jz = 1250kms"\ or = 0.0042. The average lu- 
minosity and angular diameter distances are related to 
D{z) by D^ = (l^ z)D = (1 + zfD^. 

The radial shell width, a^, is analogous to the angu- 
lar smoothing width, a^, of Section [3l The minimum 
shell width possible is set by the largest bound struc- 
tures that exist, since a regression ([2]) can only be cal- 
culated on scales over which space is expanding and a 
Hubble law is defined. Thus ccfz/H^ should be larger 
than the diameter of the largest rich clusters of galax- 
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ies, which justifies the choice made in Section [2l There 
win be similar restrictions on the choice of smoothing an- 
gle, (Jq. However, these restrictions will also be related 
to the method of angular averaging that is used, and it 
is certainly possible that angular smoothing with radial 
shells might necessit ate changes to the Gaus sian window 
averaging technique (jMcClure fc Dverll2QQ7[ ) adopted in 
Section [3l 

We will leave the development of an appropriate an- 
gular averaging scheme for future work. However, we 
note that for each shell redshift, z^, the angular correc- 
tions H{zs^O,(j)) — H{zs) will lead to corrections to the 
mean comoving distance (|3Q]) which might be expanded 
as multipoles. Equation (|3Q|) defines the monopole con- 
tribution to the distance-redshift relation. Convergence 
of the Hubble flow variance in a large data set would 
then be obtained if the dipole anisotropy converges to 
fixed value for z > Zconv^ where our preliminary inves- 
tigations suggest that the convergence scale is at least 
of order ^conv — 0.022. To consistently account for the 
residual CMB temperature dipole the residual anisotropy 
in D{z) would be up to the order of ±0.35 /i~^Mpc, with 
the exact value depending on the cosmological model. 

One intriguing question is whether an echo of the BAO 
scale could show up in a future survey on < 150 h~^Mpc 
scales if the sample statistics were greatly improved. In 
Figure [2] there are some hints of a deviation from unifor- 
mity in shells around r 100/i~^Mpc, and as remarked 
in Section 13.31 the dipole fit shows a hint of a deviation 
in the shell with mean distance f = 112 /i~^Mpc. How- 
ever, the overall effect on the average expansion is less 
than 2a in the LG frame. Since the BAO is a statisti- 
cal phenomenon, there is also a question of how it might 
manifest itself in Hubble fiow variance in our nearby vol- 
ume. If the BAO enhances wall to wall correlations, then 
we would have to account for the fact that our location 
is some distance from the closest thick wall. This might 
result in a splitting of a BAO enhancement into two fea- 
tures at Tbao i ^off, where rbao is the effective comoving 
BAO scale and Toff is the distance by which we are offset 
from the centre of the nearest wall. 

A further question which might be resolved by future 
surveys is the splitting between what in the standard cos- 
mology are known as the linear and nonlinear regimes. 
In particular, the BAO enhancement is assumed to be in 
the linear regime of perturbation theory about a FLRW 
model in the standard model. If we were to find conver- 
gence of Hubble fiow variance by ^conv — 0.022, apart 
from very small features at rbao =t ^off, then it would be 
consistent with the notion that scales z < Zconv are in the 
"nonlinear regime" while the BAO scale is in the "linear 
regime"E3- 

6.2. The minimum Hubble variance rest frame 

The multipole expansion of the Hubble fiow variance 
proposed in Section 16.11 should ideally be performed in 
the rest frame in which the radial variance in the Hubble 
flow with respect to the asymptotic global average is 

We note that such an interpretation just rehes on there ex- 
isting a scale of statistical homogeneity above which an average 
cosmological evolution can be described. It is not necessary for 
the average evolution to be described exactly by a homogeneous 
isotropic FLRW model. 



minimized. We have shown that this frame is closer to 
the LG rest frame than the CMB rest frame. The flow in 
LS frame is very close to that of the LG frame, but very 
slightly more variable. However, we have only examined 
a b oost to the LS fram e in the mean direction prescribed 
by iT ullv et aTl ()2008l ) and have not examined whether 
a slight change of the boost direction might change the 
result. 

We should also consider the possibility that the LG has 
an additional peculiar velocity with respect to the frame 
in which variance in the Hubble flow is minimized. Such 
a task is very feasible, even if computationally intensive, 
and we will report on this in future work. 

A further comparison is to independently determine 
the frame with the greatest anticorrelation between the 
residual CMB temperature dipole and the Hubble flow 
variance dipole. Does such a frame agree with the mini- 
mum Hubble variance frame, within uncertainties? 

6.3. The Hubble bubble and type la supernova 
systematics 

Type la supernovae (Snela) provide the standardiz- 
able candles which are the cornerstone of many current 
cosmological tests. The use of Snela is currently lim- 
ited by systematic uncertainties, and differences in cos- 
mological parameter estimations can be obtained when 
different light curve redu ction m ethods are used. In the 
SALT/SALT-H method (jCu^t al. 2005, 2007) empiri- 
cal light curve parameters are marginalized together with 
cosmological parameters o yer the whole da ta set, whereas 
in the MLCS2k2 method (jJha et al.ll2"007[ ) template light 
curves are determined by minimizing the distance mod- 
ulus residuals of a training set of nearby Snela, which lie 
within the range in which the Hubble flow is linear, yet 
are sufficiently distant for their peculiar velocities to be 
negligible compared with their Hubble- ffow cz. 

If the cosmic rest frame is taken to be that of mini- 
mum Hubble ffow variance on < 100 h~^Mpc scales, and 
if such a frame is close to that of the LG as our re- 
sults suggest, then an interesting systematic issue arises 
in the calibration of Snela light curves. In particular, in 
both methods one seeks to minimize the distance mod- 
ulus residuals with a nearby global linear Hubble law, 
and the standard practice is to perform such a mini- 
mization in the rest frame of the C MB rather th an th e 
rest frame of t he LG. The Union (|Kowalski et al.ll2008f ). 
Constitution (^Hicken et aTl 12009') and Union2 samples 



(jAmanullah et al.ll2010f ) contain a signiffcant number of 
data points in the rang 

43 0.015 <z< 0.02 which is be- 
low the scale ^^conv but is still conventionally deemed to 
be "within the Hubble ffow" . 

Interestingly, the redshift range 0.012 < z < 0.02 cor- 
responds to the range 36 < r < 60 /i~^Mpc over 
which the boost to the CMB rest frame was found to 
produce a smaller deviation from a uniform Hubble ffow 
than in the LG frame (c.f. Figure [2]), even though the 
Hubble ffow is signiffcantly more uniform in the LG frame 
overall. In other words, the fact that the boost to the 



20 By contrast IKessler et all (|2009t). for their full MLCS2k2 
Nearby+SDSS+SNLS+ESSENCE+HST sample, took a minimum 
redshift of z = 0.0218. There are differences in cosmological param- 
eters estimated from the Kessle r et al.l (|2009l) sample and Union, 
Constitution and Union2 samples. 
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CMB rest frame appears to best compensate for struc- 
tures in the range 30 < r < 62 h~^Mpc may have led 
to a misidentification of the minimum redshift, 2:conv, at 
which a single global linear Hubble law can be safely as- 
sumed. 

Figure [2] and Table [T] indicate that in the LG frame con- 
vergence to an almost uniform Hubble flow is achieved 
by f 65 h~^Mpc or Zconv — 0.022. This scale coin- 
cides roughly with the cutoff scale o f the Hubble bubble 
identifled in the supernovae data by Zehavi et al.' (1998) 
at z = 0.24, and conflrmed by Jha et al. (2007), using a 
MLCS2k2 sample with a reddening parameter Ry = 3.1. 
We note that over the range 60 < r < 70/i~^Mpc the 
Hubble flow is somewhat closer to uniform in the LG 
frame as opposed to the CMB frame, and the variance in 
these shel ls in either fr ame is less than th e 6.5 ± 2.2% 
found by IZehavi et al.l ([19 98). However, I Zehavi et al.l 
(|1998) and iJha et af] (^2007) worked with a far simpler 
model of Hubble flow variance in which the sample was 
divided into inner and outer sp heres. 

We can directly compare the IZehavi et al.l ()1998f ) re- 
sult with a similar analysis for the COMPOSITE sam- 
ple, by fltting a simple linear Hubble law for the 2222 
data p oint s in the interval 30 < r < 70 h~^Mpc cho- 
sen bv IZehavi et al.l (jl998f ) for their inner shell. For this 
sheh we flnd = (104.5 =b 0.6) /ikms"^ Mpc~^ in the 
CMB frame, and = (105.1 ± 0.6) /ikms-^ Mpc~^ in 
the LG frame. These values are respectively 4.40 ±0.08% 
and 4.06 ± 0.07% larger than the global asymptotic val- 
ues of determined in Section [2l They are somewhat 
lower but consisten t with the 6.5 ± 2.2% effect found by 
IZehavi et al.l ([19981 ). 

The 30 < r < 70/i~^Mpc range chosen by 
IZehavi et al.' ('1998') is equally divided between regimes 
in which the LG frame Hubble flow is closer to uniform 
in the COMPOSITE data, and alternatively in which the 
CMB frame is closer to uniform, as seen in Figure [2l This 
explains why the average values of in this range are 
closer to each other than those determined in Section [2] 
by fltting a simple linear Hubble law to the whole sample. 
The latter values, which amounted to 8.8 ± 0.2% in the 
CMB frame and 3.37 ±0.07% in the LG frame, might be 
taken as a sharper estimate of the Hubble bubble effect. 

The existence of a Hubble bubble has been contro- 
versial since as far as Snela data analysis is concerned 
the presence of the effect is dependent on the details 
of the treatnient o f extinction and reddening by dust 
(jConlev et al.ll2QQ7[ ). A Hubble bubble is found if dust in 
other galaxies has the same reddening properties as dust 
in the Milky Way but not if the re ddening parameter is 
signiflcantly reduced. iHicken e\~al. (2009) flnd no evi- 
dence for a Hubble bubble at 2; = 0.024 if the reddening 
parameter is set to Ry = 1.7. 

Our results suggest that the combination of the boost 
to the rest frame of the CMB compensating for structures 
in the range 30 < r < 62 /i~^Mpc, together with the 
treatment of parameters such as Ry as adjustable in the 
light curve reduction, may contribute signiflcantly to the 
systematic uncertainties associated with Snela. Redden- 
ing by dust in other galaxies is after all a physical quan- 
tity which should be determined independently of Snela. 
Ideally it should not be treated as a parameter which 



one can freely adjust to obtain the best flt of Hubble 
residuals. 

This issue h a s been stu died independently by 
iFinkelman et al.l (I2008L I2011D who investigated dust 
lanes in 15 E/SO galaxies and determined extinction 
properties by fltting model galaxies to the unextin- 
guished parts of the images in each of six spectral 
bands, and then subtracting these from the actual im- 
ages. They found an average val ue R^^ = 2.82 ± .38 for 
8 galaxies in their flrst study (Finkelman et al.l l2008f ). 
and Ry = 2.71 ± 0.43 for 7 galaxies in their second in- 
vestigation (jFinkelman et al.]l2011f ). For the combined 
sample Ry = 2.77 ± 0.41. This value is a little lower 
than the Milky Way value Ry = 3.1 but consistent with 
it within the uncertainty. 

Our results suggest that the convergence scale ^conv — 
0.022 is clos e to that of t he H ubble bubble originally 
proposed by IZehavi et al.l (|1998l ) , but the magnitude of 
the Hubble bubble effect is smaller when viewed in the 
LG frame. For consistency MLCS2k2 Snela data should 
be reduced using Ry values consistent with independent 
determinati ons, e.g., R^^ = 2 .77 ± 0.41 as suggested by 
the work of IFinkelman et all (|2008l . I2011f ). As discussed 
bv ISmale Wiltshird ([20Tlh this is also important for 
cosmological model comparison. 

Since the First Amendment Snela data of 
iTurnbull et all (|20T2h was reduced with Ry = 1.7, 
it would be interesting to know to what extent the 
difference in the amplitude of the bul k flow velocity 
from that of the COMPOSITE sample of (|Watkins et al.l 
I2OO9I ) is due to the choic e of the R^^ pa rameter. We 
suggest that the data set of ITurnbull et al.l (|2012l ) should 
be reanalysed with Ry = 2.77 ±0.41, and by the method 
of Section [2] in the LG frame. 



6.4. The asymptotic global Hubble constant 

The COMPOSITE sample enables us to determine the 
relative Hubble flow, but does not constrain the overall 
normalization of the distance scale and consequently the 
precise value of the global asymptotic Hubble constant. 

The Hubble constant has recently been determined to 
high accuracy by the SHOES survey as = 73.8 ± 

2.4kms~^ Mpc~^ (jRiess et al.l l201 ih . Independent es- 
timates of the Hubble constant us ing; baryon acoustic 
oscillations at a variety of redshifts (iPercival et al.l 120101 : 
iBeutler et all 120111 : iBlake et al.l 1201 ID . have vielded val- 
ues = {68.2 ±2.2, 67 ±3.2, 68.1 ± 1.7} kms-^Mpc"^ 
respectively in the (possibly curved) AC DM model. 
While these estimates are consistent with the lRiess et al.l 
(|2011l ) result at the 2a level, a further increase in pre- 
cision could lead to tension. A signiflcant difference be- 
tween the baryonic acoustic oscillation method and the 
SHOES survey is that the former relies on the flt to the 
FLRW model at larger redshifts 0.1 < z < 0.6, whereas 
the latter is less model dependent but relies on establish- 
ing a cosmic distance ladder of distance indicators such 
as Cepheid variables on very nearby scales. 

If we identify the cosmic rest frame with that of mini- 
mum Hubble flow variance, then the impact of perform- 
ing all cosmological tests in such a frame rather than in 
the CMB frame needs to be carefully considered. The im- 
pact is likely to be most signiflcant on those tests which 
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directly use measurements on z < 0.022 scales. Whether 
this has an impact on measurements that establish the 
cosmic distance ladder is an intriguing question which 
should be investigated once the minimum Hubble flow 
variance rest frame is positively identified. 

6.5. Large angle CMB anomalies 

There are several observations concerning the large an- 
gle multipoles of the CMB anisotropy spectrum, which 
may be considered anomalous to varying degrees of 
statistical significance. These include: (i) the power 
asymme try between the northern and southern hemi- 
spheres (|Tegm ark et al.l l2003l : lEnksen et aT 2004, 2007) ; 
(ii ) the low q u adrupole po wer (Tegmark et al. 2003'; 
Ide Oliveira- Costa et al.ll2004f): fiii) the alignment of the 
quadrupole and octupole (|de O liveira-Costa et al. 2004; 
Schwarz et all 120041 : iLand fc M agueiio 2005; Copi et al. 
2006"); and (iv) the parity asymmetry (^Kim fc Naselskvt 
2010). It is beyond the scope of the present paper to 
investigate all these anomalies. However, it is clear that 
our proposal to revisit a significant feature of the CMB 
anisotropy analysis, namely the nature of the dipole, will 
introduce systematics which would necessitate a reassess- 
ment of all of these issues. 

There are two obvious potential lines of inquiry: 

• The propagation of photons through the fore- 
grounds contributing to the Hubble fiow variance 
may produce a multipole signature which differs 
subtly from the pure dipole signature ([2T]) associ- 
ated with a Lorentz boost; 

• Since the dipole subtraction is an integral part of 
the map-making procedure, differences in dipole 
subtraction may lead to subtle differences i n the 
cleaning of galactic foregrounds (jTegmark et al.l 
[20031 

A study by iFreeman et al.l (|2006[ ) found that of several 
possible systematic errors, a 1-2% error in the CMB 
dipole subtraction stood out as being an effect which 
could potentially resolve the power asymmetry anomaly. 

We note that while a 1-2% change in the dipole would 
not affect the power on small angles, its effect on the 
large angle multipoles would require a redrawing of the 
CMB sky maps. Any such redrawing could also alter 
other large angle features, such as the Cold Spot. 

6.6. The dark flow 

The determination of peculiar v elocities via the 
kinematic Sunyaev-Zel'dovich effect fKashlinskv et all 
^009a,b, 2010; ^Hand et al. 2011; |Lavaux et al. 20121) is 
particularly interesting, since it is a method which uses 
the CMB anisotropy spectrum, rather than directly 
looking at the distance-redshift relation. Furthermore, 
the very large bulk fiow tha t has been claimed by 
iKashlinskv et al.l ()2009ai l2010f ) is controversial and has 
not be e n reproduced in subs e quent studie s , e.g., iKeislerl 
(|2009h : ICMTorne et al.l (|20TTh : iHand eA~^ m m. 

A significant featur e of the measurements of 
IKashlinskv et al.l (j2009al lbl. I2QToI ) is the claim that the re- 
sult is independent of many systematics, since both fore- 
ground and cosmological dipoles and quadrupoles have 
been subtracted in a consistent way. It is claimed that 



the peculiar velocities inferred are those of the galaxy 
clusters with respect to the CMB in their own rest 
frames. Nonetheless, since the direction of the reported 
bulk "dark flow", (^,6) = (267°, 234°) ± 15° coincides 
so closely with that of the residual CMB dipole in the 
LG rest frame, all potential systematic errors must be 
carefully considered. 

Although we have not gone throu gh all of the de- 
tails of the highly technical analysis of IKashlinskv et aTl 
(|2009b l we can identify one potential systematic error 
in the reported description of their procedure. To de- 
termine the kinem atic Sunyaev-Zel'dovich (kSZ) dipole 
IKashlinskv et al.l (|2009 b) must isolate it from other con- 
tributions in their equation (1). In particular, they sub- 
tract the contribution of the thermal Sunyaev-Zel'dovich 
effect. The calculation of the thermal Sunyaev-Zel'dovich 
contribution involves a determination of the temperature 
of the thermal plasma in the galaxy cluster. This is de- 
termined iteratively using empirical astrophysical rela- 
tions together with the cluster redshift. Total rest-frame 
[0.1-2.4] keV band cluster luminosities are then deter- 
mined from the recalculated fluxes using a standard con- 
version with the cosmological luminosity distance and a 
tempe rature-dependent K-correction (Kas hlinskv et al.l 
l2009bf ). All of these calculations have been performed 
using a perfectly isotropic distance-redshift relation, as 
one expects in the standard cosmology. 

If the actual distance-redshift relation has a small 
dipole anisotropy due to foreground inhomogeneities 
then it is possible that a small systematic offset with 
a dipole dependence is introduced in the cluster temper- 
ature estimate. Altern atively, we i night s ay given their 
methodology the Kashl inskv et al.l ()2009a[ ) result is con- 
sistent with our hypothesis. However, it is due to a differ- 
ential expansion of space on < 65 /i~^Mpc scales, rather 
than a bulk flow extrapolated to the scale of the sample 
measurement. This should be further investigated. 

If we are correct, then investigations of the kSZ effect 
which employ a d ifferent methodology, such as that of 
I Hand et al.l ()2012f ) based on correlations between cluster 
pa irs, will not see bulk flows of t he magnitude claimed 
bv IKashlinskv et al.l ([20095 l20Tol ). The kSZ effect is of 
course due to the local CMB dipole at the cluster loca- 
tion. With our interpretation, this temperature dipole 
will include contributions from both a peculiar velocity 
and from the differential expansion of space due to in- 
homogeneities in the vicinity of the galaxy cluster. Es- 
timates of the maximum Hubble flow variance based on 
void/wall statistics should therefore put bounds on the 
magnitude of what is purely assigned to a "peculiar ve- 
locity" in the standard framework. 

6.7. Conclusion 

While much work remains to be done, our results sug- 
gest that a fundamental revision of the treatment of 
peculiar velocities may shed new light on many of the 
puzzles raised by bulk flows, and perhaps even resolve 
some of the associated anomalies. If a large fraction of 
the CMB dipole is due to a residual anisotropy in the 
distance-redshift relation, as our results seem to suggest, 
then this will also have important consequences for un- 
derstanding large-angle features in the CMB anisotropy 
spectrum. 

Peculiar velocities will always play a role in observa- 
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tional cosmology - galaxies in clusters exhibit peculiar 
velocities with respect to the cluster barycentres, and 
this is directly observed in the "fingers of god" redshift 
space distortions. However, the most natural conclusion 
of our work is that on scales larger than gravitationally 
bound systems the variance of the Hubble flow should be 
treated as the differential expansion of regions of vary- 
ing density, which have decelerated by different amounts 
from the initial uniform distribution at the epoch of last 
scattering. One should only talk about "infall" if the 
physical distance between objects is actually decreasing 
with time, rather than applying it to the case of denser 
regions which are expanding less than the average. While 
directional forces are the basis of Newtonian mechanics 
in Euclidean space, there is nothing in general relativity 
which demands that such notions should apply to scales 
of tens of megaparsecs over which space is expanding. 

In recent years there has been some discussion about 
whether it is conceptually more correct to think of space 
as expanding, or whether the treatment of the expan- 
sion by a simple Doppler law on a fixed backg:round is 
sufficient (|Abramowicz et al.l I2QQ7I : IChodorowskil I2QQ7I : 



IGr0n & Ele:ar0y| 120071; "Francis et al.l 120071; iLewis elalj 
2007; Bunn & Hose 2009; Roukema 2009; Faraoni 2010|). 
In particular, lAbramowicz et al.. (^07) showed that the 
expansion of space can in principle be observatio nally 
determined. Some other authors, e.g., IBunn fc HoggI 
(2009), maintained that the Doppler law picture is still 
useful. The debate involved thought experiments con- 
ducted within homogeneous isotropic cosmological mod- 
els. Our results suggest that, as far as actual observations 
are concerned, variance in the Hubble law over scales of 
tens of megaparsecs cannot be simply reduced to a boost 
at a point; space really is expanding, and by differential 
amounts. 
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APPENDIX 
DIPOLE ESTIMATION 
2 



In Section 13.21 for each radial shell we minimize = Xli ^^j^ i^i ~ — 1^ ■ ni) , where /3 = {Px^ Py^ Pz) = 

{P cos id cos bd^ /3cos£(i sin6(i, P sin bd)^ = {uxi, Uyi^ n^i) = (cos -^^ cos 6^, cos sin 6^, sin^^). Hi = czi/ri and its 
uncertainty is given by (p!3|) . Minimization with respect to the four independent parameters = {Hd, Px^ Py, Pz) 
yields the linear system 



Hd 
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^zi ^xi 



Hi n 



yi 



H, 



Hi 



(J 



H, 



Hi 



i H, 
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Hi 



which is readily solved. The covariance matrix for the original variables = {H^, /3,£^,bd) is obtained straightfor- 
wardly from the covariance matrix Cov(X'^, X'^) by the standard relation 



Cov(r°,r'') = ^Cov(x^x'') 



dXd 



(A5) 



where dY^ /dX^ is the relevant Jacobian. 



REFERENCES 



Abramowicz, M. A., Bajtlik, S., Lasota, J. P. &; Moudens, A., 

2007 Acta Astron. 57, 139. 
Alnes, H. & Amarzguioui, M., 2006 Phys. Rev. D 74, 103520. 
Amanullah, R., Lidman, C, Rubin, D., et al., 2010 ApJ 716, 712. 
Antoniou, I. &; Perivolaropoulos, L., 2010 J. Cosmol. Astropart. 

Phys., JCAP12(2010)012. 
Bennett, C. L., Halpern, M., Hinshaw, G., et a/., 2003, ApJS 148, 

1. 

Rentier, P., Blake, C, Colless, M., et al, 2011 MNRAS 416, 3017. 
Bilicki, M., Chodorowski, M., Mamon, G. A. & Jarrett, T., 2011 
ApJ 741, 31. 

Blake, C., Kazin, E., Beutler, P., et al, 2011 MNRAS 418, 1707. 
Buchert, T., 2000 Gen. Relativ. Grav. 32, 105. 



Buchert, T., 2008 Gen. Relativ. Grav. 40, 467. 

Bunn, E. F. & Hogg, D. W., 2009 Am. J. Phys. 77, 688. 

Chodorowski, M., 2007 MNRAS 378, 239. 

Cohn, J., Mohayaee, R., Sarkar, S. & Shafieloo, A., 2011 MNRAS 
414, 264. 

Conley, A., Carlberg, R. G., Guy, J., Howell, D. A., Jha, S., 

Riess, A. G. &; Sullivan, M., 2007 ApJ 664, L13. 
Copi, C. J., Huterer, D., Schwarz, D. J. &; Starkman, G. D., 2006 

MNRAS 367, 79. 
Dai, D. C, Kinney, W. H. & Stojkovic, D., 2011 J. Cosmol. 

Astropart. Phys., JCAP04(2011)015. 
de Oliveira-Costa, A., Tegmark, M., Zaldarriaga, M. &; Hamilton, 

A., 2004 Phys. Rev. D 69, 063516. 



20 



Wiltshire, Smale, Mattsson & Watkins 



Erdogdu, P., Lahav, O., Huchra, J. P., et a/., 2006 MNRAS 373, 
45. 

Eriksen, H. K., Hansen, F. K., Banday, A. J., Gorski, K. M. &; 

Lilje, P. B., 2004 ApJ 605, 14; (E) 609, 1198. 
Eriksen, H. K., Banday, A. J., Gorski, K. M., Hansen, F. K. & 

Lilje, P. B., 2007 ApJ 660, L81. 
Faraoni, V., 2010 Gen. Relativ. Grav. 42, 851. 
Feldman, H. A., Watkins, R. & Hudson, M. J., 2010 MNRAS 

407, 2328. 

Finkelman, I., Brosch, N., Kniazev, A. Y., et al., 2008 MNRAS 
390, 969. 

Finkelman, I., Brosch, N., Kniazev, A. Y., et al, 2011 MNRAS 
409 727. 

Fixsen, D. J., Cheng, E. S., Gales, J. M., Mather, J. C., Shafer, 

R. A. & Wright, E. L., 1996, ApJ 473, 576. 
Francis, M. J., Barnes, L. A., James, J. B. & Lewis, G. F., 2007 

PASA 24, 95. 

Freedman, W. L., Madore, B. F., Gibson, B. K., et al, 2001 ApJ 
553, 47. 

Freeman, P. E., Genovese, C. R., Miller, C. J., Nichol R. C. & 

Wasserman, L., 2006 ApJ 638, 1. 
Gordon, C., Land, K. &; Slosar, A., 2008 MNRAS 387, 371. 
Gr0n, 0. &; Elgar0y, 0., 2007 Am. J. Phys. 75, 151. 
Guy, J., Astier, P., Nobih, S., Regnault, N. & Pain, R., 2005 

A&A 443, 781. 
Guy, J., Astier, P., Baumont, S., et a/., 2007 A&A 466, 11. 
Hand, N., Addison, G. E., Aubourg, E., et al, 2012 Phys. Rev. 

Lett. 109, 041101. 
Hicken, M., Wood-Vasey, W. M., Blondin, S., et al, 2009 ApJ 

700, 1097. 

Hogg, D. W., Eisenstein, D. J., Blanton, M. R., Bahcall, N. A., 
Brinkmann, J., Gunn, J. E. &; Schneider, D. P., 2005 ApJ 624, 
54. 

Hoyle, F. & Vogeley, M. S., 2002 ApJ 566, 641; 
Hoyle, F. & Vogeley, M. S., 2004 ApJ 607, 751. 
Hudson, M. J., Smith, R. J., Lucey, J. R. & Branchini, E., 2004 

MNRAS 352, 61. 
Iwata, I. & Chamaraux, P., 2011 A&A 531, A87. 
Jha, S., Riess, A. G. & Kirshner, R. P., 2007 ApJ 659, 122. 
Karachentsev, I. D., Kashibadze, O. G., Makarov, D. I. &: Tully, 

R. B., 2009 MNRAS 393, 1265. 
Kashlinsky, A., Atrio-Barandela, F., Kocevski, D. &; Ebeling, H., 

2009a ApJ 686, L49. 
Kashlinsky, A., Atrio-Barandela, F., Kocevski, D. & Ebeling, H., 

2009b ApJ 691, 1479. 
Kashlinsky, A., Atrio-Barandela, F., Ebeling, H., Edge, A. &; 

Kocevski, D., 2010 ApJ 712, L81. 
Kaiser, N., 1987 MNRAS 227, 1. 

Kass, R.E. & Raftery, A.E., 1995 J. Am. Statist. Assoc. 90, 773. 
Keisler, R., 2009 ApJ 707, L42. 

Keisler, R., Reichardt, C. L., Aird, K. A., et al, 2011 ApJ 743, 
28. 

Kessler R., Becker, A., Cinabro, D., et al, 2009 ApJS 185, 32. 

Kim, J. & Naselsky, P., 2010 ApJ 714, L265. 

Kocevski, D. D. & Ebehng, H., 2006 ApJ 645, 1043. 

Kogut, A., Lineweaver, C., Smoot, G. F., et al, 1993 ApJ 419, 1. 

Kowalski, M., Rubin, D., Aldering, G., et a/., 2008 ApJ 686, 749. 

Land, K. & Magueijo, J., 2005 Phys. Rev. Lett. 95, 071301 



Lauer, T. R. &; Postman, M., 1994, ApJ 425, 418. 
Lavaux, G., Tully, R. B., Mohayaee, R. &; Colombi, S., 2010 ApJ 
709, 483. 

Lavaux, G., Afshordi, N. and Hudson, M. J., 2012 
arXiv: 1207. 1721 

Lewis, G. F., Francis, M. J., Barnes, L. A. &; James, J. B., 2007 

MNRAS 381, L50. 
Li, N. & Schwarz, D. J., 2008 Phys. Rev. D 78, 083531. 
Linder, E. V., 2005 Phys. Rev. D 72, 043529. 
McClure, M. L. & Dyer, C. C., 2007 New Astron. 12, 533. 
Osborne, S. J., Mak, D. S. Y., Church, S. E. & Pierpaoh, E., 2011 

ApJ 737, 98. 

Pan, D. C, Vogeley, M. S., Hoyle, F., Choi, Y. Y., & Park, C, 

2012 MNRAS 421, 926. 
Partridge, R. B. & Wilkinson, D. T., 1967 Phys. Rev. Lett. 18, 

557. 

Peebles, P. J. E., 1993 Principles of Physical Cosmology, 

(Princeton University Press). 
Peebles, P. J. E. &; Wilkinson, D. T., 1968 Phys. Rev. 174, 2168. 
Percival, W., Reid, B. A., Eisenstein, D. J., et a/., 2010 MNRAS 

401, 2148. 

Press, W. H., Flannery, B. P., Teukolsky, S. A. & Vetterling, 
W. T., 1986 Numerical Recipes, (Cambridge University Press). 

Riess, A. G., Macri, L., Casertano, S., et al, 2011 ApJ 730, 119; 
(E) 732, 129. 

Roukema, B. F., 2010 MNRAS 404, 318. 

Sandage, A., 1986 ApJ 307, 1. 

Schwarz, D. J., Starkman, G. D., Huterer, D. and Copi, C. J., 

2004 Phys. Rev. Lett. 93, 221301. 
Smale, P. R., 2012 PhD thesis. University of Canterbury. 
Smale, P. R. & Wiltshire, D. L., 2011 MNRAS 413, 367. 
Springob, C. M., Masters, K. L., Haynes, M. P., et al, 2007 ApJS 

172, 599. 

Sylos Labini, F., Vasilyev, N. L., Pietronero, L. &; Baryshev 

Y. v., 2009 Europhys. Lett. 86, 49001. 
Tegmark, M., de Oliveira-Costa, A. &; Hamilton, A., 2003 Phys. 

Rev. D 68, 123523. 
Thomas, S. A., Abdalla, F. B. &; Lahav, O., 2011 Phys. Rev. 

Lett. 106, 241301. 
Trotta, R., 2007 MNRAS 378, 72. 

Tully, R. B., Shaya, E. J., Karachentsev, L D., Courtois, H. M., 
Kocevski, D. D., Rizzi, L. & Peel, A., 2008 ApJ 676, 184. 

Turnbull, S. J., Hudson, M. J., Feldman, H. A., Hicken, M., 
Kirshner, R. P. &; Watkins, R., 2012 MNRAS 420, 447. 

Watkins, R., Feldman, H. A. & Hudson, M. J., 2009 MNRAS 
392, 743. 

Weyant, A., Wood-Vasey, M., Wasserman, L. & Freeman, P., 

2011 ApJ 732, 65. 
Wiltshire, D. L., 2007a New J. Phys. 9, 377. 
Wiltshire, D. L., 2007b Phys. Rev. Lett. 99, 251101. 
Wiltshire, D. L., 2008 Phys. Rev. D 78, 084032. 
Wiltshire, D. L., 2009 Phys. Rev. D 80, 123512. 
Wiltshire, D. L., 2011 Class. Quantum Grav. 28, 164006. 
Zehavi I., Riess A. G., Kirshner R. P. & Dekel A. 1998 ApJ 503, 

483. 



Hubble flow variance and the cosmic rest frame 



21 




(c) 3: 25 - 37.5 h'^ Mpc = 514. (d) 4: 37.5 - 50 h'^ Mpc = 731. 




(g) 7: 75 - 87.5 h'^ Mpc = 414. (h) 8: 87.5 - 100 h'^ Mpc = 304. 




(i) 9: 100 - 112.5 h'^ Mpc N = 222. (j) 10: 112.5 - 156.25 h'^ Mpc A^ = 280. 







m 










*• 













(k) 11: 156.25 - 417.4 h'^ Mpc A^ = 91. 



Fig. 1. — The angular distribution of the individual data points in the COMPOSITE sample is plotted in a manner similar to Figure 1 of 
Ipteldman et al. ( 2010), but distributed into the (unprimed) radial shells of Table [T] The radius of each data point indicates the magnitude of 
the peculiar velocity relative to HsVi, where Hs is the value given for the LG frame in Tabled in the unprimed case. Positive and negative 
peculiar velocities are colour-coded red and blue (or darker and lighter in greyscale). In all figures, the galactic longitudes ^ = 0°, 180°, 360° 
are on the right edge, centre and left edge respectively. The same scale is used in all subsequent figures. 
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TABLE 1 

Hubble flow variance in radial shells in CMB and LG frames, for two different choices of 

SHELLS, THE SECOND CHOICE LABELLED BY PRIMES WITH A COMMON FINAL SHELL. IN EACH CASE 

rs < r < Ts+i. We tabulate the inner shell radius, Ts] the weighted mean distance, fs] THE 
SHELL Hubble constants, (^s)^^^ and {Hs)i^q in the CMB and LG frames, and their 

UNCERTAINTIES DETERMINED BY LINEAR REGRESSION WITHIN EACH SHELL, TOGETHER WITH ITS "GOODNESS 

OF fit" probability Qs and reduced (for u = Ns - 1); \nB WHERE B IS THE Bayes factor for 

THE RELATIVE PROBABILITY THAT THE LG FRAME HAS MORE UNIFORM SHs = THAN THE CMB FRAME 
WHEN IS SUMMED IN ALL SHELLS WITH r > Ts- Hg AND Q-g ARE GIVEN IN UNITS /lKMS~^ MPC""*-. 



Shell s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Ns 


92 


505 


514 


731 


819 


562 


414 


304 


222 


280 


91 




2.00 


12.50 


25.00 


37.50 


50.00 


62.50 


75.00 


87.50 


100.00 


112.50 


156.25 


fs {h-^Mpc) 


5.43 


16.33 


30.18 


44.48 


55.12 


69.24 


81.06 


93.75 


105.04 


126.27 


182.59 


(-^«)CMB 


173.9 


111.1 


110.4 


104.1 


102.7 


103.8 


102.1 


102.8 


104.4 


102.1 


100.1 




6.8 


1.5 


1.1 


0.7 


0.7 


0.8 


0.9 


0.9 


1.0 


0.8 


1.7 


(Qs)cMB' 


0.000 


0.000 


0.000 


0.037 


0.985 


0.997 


1.000 


1.000 


1.000 


0.979 


0.999 


(XV^)CMB 


24.639 


5.802 


1.741 


1.096 


0.896 


0.841 


0.593 


0.604 


0.630 


0.835 


0.581 


(^^)lg 


117.9 


103.1 


106.5 


105.5 


104.8 


102.1 


102.8 


103.2 


103.7 


102.4 


101.0 




4.6 


1.4 


1.0 


0.7 


0.7 


0.7 


0.9 


0.9 


1.0 


0.8 


1.7 


(Qs)lG' 


0.000 


0.000 


0.000 


0.000 


0.998 


0.940 


1.000 


1.000 


1.000 


0.993 


0.999 


(xV-)lg 


23.656 


7.767 


2.185 


1.419 


0.864 


0.909 


0.594 


0.542 


0.622 


0.803 


0.590 


\nB {r> rs) 


58.62 


16.95 


8.43 


1.71 


1.98 


2.71 


1.64 


1.78 


1.67 


0.49 





Shell s 


r 


2' 


3' 


4' 


5' 


6' 


7' 


8' 


9' 


10' 


11 


Ns 


321 


513 


553 


893 


681 


485 


343 


273 


164 


206 


91 


rs (h-^Mpc) 


6.25 


18.75 


31.25 


43.75 


56.25 


68.75 


81.25 


93.75 


106.25 


118.75 


156.25 


fs (h-^Mpc) 


12.26 


23.46 


37.61 


49.11 


61.74 


73.92 


87.15 


99.12 


111.95 


131.49 


182.59 




120.8 


108.8 


105.9 


103.6 


104.2 


102.6 


102.7 


103.8 


102.0 


102.3 


100.1 


(^^)CMB 


2.1 


1.2 


0.9 


0.7 


0.8 


0.8 


0.9 


0.9 


1.0 


0.9 


1.7 




0.000 


0.000 


0.000 


0.555 


0.959 


1.000 


1.000 


1.000 


0.992 


0.997 


0.999 


(XV^)CMB 


9.496 


2.506 


1.421 


0.993 


0.908 


0.633 


0.624 


0.658 


0.754 


0.745 


0.581 




103.5 


103.5 


103.9 


106.6 


103.9 


102.0 


103.2 


103.6 


101.6 


102.7 


101.0 


^^)lg 


1.8 


1.1 


0.9 


0.7 


0.8 


0.8 


0.9 


0.9 


1.0 


0.9 


1.7 




0.000 


0.000 


0.000 


0.031 


0.960 


1.000 


1.000 


1.000 


0.996 


0.999 


0.999 


(xV-)lg 


11.427 


3.246 


1.792 


1.090 


0.907 


0.701 


0.592 


0.608 


0.728 


0.711 


0.590 


\nB {r> rs) 


30.09 


8.99 


3.22 


0.92 


2.47 


1.68 


1.37 


1.30 


0.64 


0.39 





TABLE 2 

The analysis of Table [T] is repeated removing all points which contribute a value x^ > 5 in 
both the cmb and lg frames in both the primed and unprimed shells. a total 4398 points 

REMAIN. 



Sheh s 12345678 9 10 11 

Ns 78 457 494 713 799 555 413 303 221 274 91 

rs {h-^Mpc) 2.00 12.50 25.00 37.50 50.00 62.50 75.00 87.50 100.00 112.50 156.25 

fs (h-^Mpc) 7.23 16.22 30.12 44.52 55.11 69.25 81.06 93.68 105.10 126.51 182.59 

(i/s)cMB 1^2.7 109.5 108.2 103.5 101.4 103.0 102.0 103.1 104.1 102.1 100.1 

(^«)CMB 1-0 

Ws)cMB' 0-000 0.000 0.024 1.000 1.000 1.000 1.000 1.000 1.000 1.000 0.999 

(X^/^)CMB 9-092 1.818 1.130 0.739 0.606 0.664 0.580 0.577 0.603 0.685 0.581 

(Hs)i^Q 107.8 98.2 103.7 105.4 103.6 101.4 102.7 103.5 103.4 102.4 101.0 

(as)LG 1-0 0.8 0.7 0.7 0.9 0.9 1.0 0.8 1.7 

(Q3)lG' 0.000 0.000 0.000 0.980 1.000 1.000 1.000 1.000 1.000 1.000 0.999 

(x^/^)lG 4.048 2.340 1.431 0.894 0.583 0.723 0.578 0.513 0.595 0.667 0.590 

\nB(r>rs) 50.84 11.42 5.97 1.04 1.58 2.16 1.53 1.67 1.52 0.44 

Sheh s V 2' 3' 4' 5' 6' 7' 8' 9' 10' 11 

Ns 

rs (h-^Mpc) 
fs (h-^Mpc) 

(:^«)cMB 

(^«)CMB 
(Qs)cMB' 
(XV^)CMB 

(^«)lg 
(Qs)lg^ 
(xV-)lg 

InB {r> rs) 



284 


488 


532 


869 


669 


481 


343 


271 


160 


204 


91 


6.25 


18.75 


31.25 


43.75 


56.25 


68.75 


81.25 


93.75 


106.25 


118.75 


156.25 


12.25 


24.05 


37.47 


49.17 


61.75 


73.92 


87.15 


99.13 


111.81 


131.52 


182.59 


119.5 


107.3 


105.3 


102.4 


103.0 


102.2 


102.7 


103.9 


101.7 


102.3 


100.1 


2.2 


1.2 


1.0 


0.7 


0.8 


0.8 


0.9 


0.9 


1.1 


0.9 


1.7 


0.000 


0.000 


0.649 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.999 


4.042 


1.313 


0.975 


0.645 


0.638 


0.587 


0.624 


0.604 


0.561 


0.690 


0.581 


99.7 


101.1 


103.6 


105.7 


102.7 


101.6 


103.2 


103.6 


101.5 


102.7 


101.0 


1.8 


1.1 


0.9 


0.7 


0.8 


0.8 


0.9 


0.9 


1.0 


0.9 


1.7 


0.000 


0.000 


0.001 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.999 


3.581 


1.496 


1.197 


0.699 


0.653 


0.650 


0.592 


0.552 


0.563 


0.655 


0.590 


23.97 


6.29 


2.31 


0.59 


1.86 


1.53 


1.32 


1.24 


0.57 


0.42 
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TABLE 3 

The analysis of Table [T] is repeated removing all points which contribute a value > 5 in 
either the cmb or lg frames in both the primed and unprimed shells. a total 4212 points 

REMAIN. 



Shell s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Ns 


56 


385 


454 


692 


788 


548 


410 


301 


217 


270 


91 


rs (h-^Mpc) 


2.00 


12.50 


25.00 


37.50 


50.00 


62.50 


75.00 


87.50 


100.00 


112.50 


156.25 


fs {h-^Mpc) 


9.04 


17.87 


30.33 


44.58 


55.10 


68.74 


80.95 


93.74 


105.08 


126.20 


182.59 




155.4 


109.3 


107.6 


103.0 


101.4 


103.9 


102.0 


103.4 


104.6 


102.6 


100.1 




5.0 


1.6 


1.2 


0.8 


0.7 


0.8 


0.9 


0.9 


1.0 


0.8 


1.7 


(Qs)cMB' 


0.000 


0.228 


0.991 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.999 


(X /^)CMB 


2.895 


1.053 


0.850 


0.632 


0.546 


0.616 


0.550 


0.538 


0.498 


0.613 


0.581 


(^«)lg 


114.3 


101.3 


104.7 


105.0 


103.5 


102.8 


102.9 


103.6 


103.7 


102.8 


101.0 




3.6 


1.5 


1.1 


0.8 


0.7 


0.8 


0.9 


0.9 


1.0 


0.8 


1.7 




0.000 


0.529 


0.778 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.999 


(xV-)lg 


2.275 


0.993 


0.949 


0.745 


0.523 


0.647 


0.549 


0.482 


0.519 


0.606 


0.590 


InB (r > rs) 


50.81 


10.83 


5.31 


1.78 


2.40 


3.02 


1.99 


2.19 


1.98 


0.59 





Shell s 

Ns 

rs (h-^Mpc) 
fs (h-^Mpc) 

(^«)CMB 

(Qs)cMB' 

(XV^)CMB 

^^)lg 
(xV-)lg 

InB {r> rs) 



r 

223 
6.25 
12.92 
124.5 

2.4 
0.000 
2.414 
105.2 

2.0 
0.000 
1.444 
29.83 



2' 
436 
18.75 
24.39 
106.4 
1.3 
0.994 
0.839 
101.3 
1.2 
0.863 
0.926 
5.66 



3' 
501 
31.25 
37.80 
105.7 
1.0 
1.000 
0.794 
104.4 
1.0 
0.972 
0.883 
2.75 



4' 
856 
43.75 
49.23 
102.3 
0.7 
1.000 
0.588 
105.3 
0.7 
1.000 
0.616 
1.08 



5' 
657 
56.25 
61.51 
102.9 
0.8 
1.000 
0.568 
102.9 
0.8 
1.000 
0.566 
2.29 



6' 
477 
68.75 
73.84 
103.0 
0.8 
1.000 
0.556 
102.9 
0.8 
1.000 
0.588 
2.01 



7' 
339 
81.25 
87.08 
103.0 
0.9 
1.000 
0.568 
103.5 
0.9 
1.000 
0.539 
1.66 



8' 
268 
93.75 
99.03 
104.0 
0.9 
1.000 
0.554 
103.6 
0.9 
1.000 
0.506 
1.59 



9' 
158 
106.25 
111.89 
101.9 
1.1 
1.000 
0.486 
101.6 
1.1 
1.000 
0.510 
0.83 



10' 
201 
118.75 
131.38 
103.1 
0.9 
1.000 
0.610 
103.4 
0.9 
1.000 
0.591 
0.62 



11 
91 
156.25 
182.59 
100.1 
1.7 
0.999 
0.581 
101.0 
1.7 
0.999 
0.590 



TABLE 4 

OVING ALL F 

EITHER THE CMB OR LG FRAMES IN EITHER THE PRIMED OR UNPRIMED SHELLS. A TOTAL 4156 POINTS 

REMAIN. 



Sheh s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Ns 


30 


366 


447 


691 


786 


547 


410 


301 


217 


270 


91 


rs (h-^Mpc) 


2.00 


12.50 


25.00 


37.50 


50.00 


62.50 


75.00 


87.50 


100.00 


112.50 


156.25 


fs {h-^Mpc) 


8.44 


18.32 


30.20 


44.58 


55.10 


68.74 


80.95 


93.74 


105.08 


126.20 


182.59 


(;^«)CMB 


138.4 


107.2 


107.1 


103.0 


101.2 


103.8 


102.0 


103.4 


104.6 


102.6 


100.1 


(^«)CMB 


6.3 


1.6 


1.2 


0.8 


0.7 


0.8 


0.9 


0.9 


1.0 


0.8 


1.7 


(Qs)cMB' 


0.320 


0.992 


0.997 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.999 


(XV^)CMB 


1.103 


0.829 


0.823 


0.627 


0.531 


0.611 


0.550 


0.538 


0.498 


0.613 


0.581 




106.9 


101.7 


104.2 


104.9 


103.4 


102.7 


102.9 


103.6 


103.7 


102.8 


101.0 


^^)lg 


4.9 


1.5 


1.1 


0.8 


0.7 


0.8 


0.9 


0.9 


1.0 


0.8 


1.7 




0.464 


0.530 


0.907 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.999 


(xV-)lg 


1.001 


0.993 


0.912 


0.738 


0.519 


0.638 


0.549 


0.482 


0.519 


0.606 


0.590 


\nB {r> rs) 


21.00 


8.00 


4.88 


1.74 


2.33 


2.96 


1.99 


2.19 


1.98 


0.59 





Sheh s 

Ns 

rs (h-^Mpc) 

fs (h-^Mpc) 

(;^s)cMB 

(^«)CMB 

(Qs)cMB' 

(XV^)CMB 

^^)lg 
(xV-)lg 

\nB {r> rs) 



r 

181 
6.25 
14.27 
112.5 

2.3 
0.504 
0.995 
103.0 

2.1 
0.243 

I. 071 

II. 34 



2' 
429 
18.75 
24.55 
106.3 
1.3 
0.999 
0.803 
100.9 
1.2 
0.966 
0.879 
4.80 



3' 
498 
31.25 
37.83 
104.9 
1.0 
1.000 
0.770 
104.3 
1.0 
0.980 
0.875 
2.07 



4' 
854 
43.75 
49.22 
102.2 
0.7 
1.000 
0.578 
105.2 
0.7 
1.000 
0.610 
1.03 



5' 
655 
56.25 
61.49 
102.7 
0.8 
1.000 
0.555 
102.7 
0.8 
1.000 
0.555 
2.22 



6' 
477 
68.75 
73.84 
103.0 
0.8 
1.000 
0.556 
102.9 
0.8 
1.000 
0.588 
2.01 



7' 
339 
81.25 
87.08 
103.0 
0.9 
1.000 
0.568 
103.5 
0.9 
1.000 
0.539 
1.66 



268 
93.75 
99.03 
104.0 

0.9 
1.000 
0.554 
103.6 

0.9 
1.000 
0.506 
1.59 



9' 
158 
106.25 
111.89 
101.9 
1.1 
1.000 
0.486 
101.6 
1.1 
1.000 
0.510 
0.83 



10' 
201 
118.75 
131.38 
103.1 
0.9 
1.000 
0.610 
103.4 
0.9 
1.000 
0.591 
0.62 



11 
91 
156.25 
182.59 
100.1 
1.7 
0.999 
0.581 
101.0 
1.7 
0.999 
0.590 
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(a) 



40 60 80 100 
(ts) (h-^ Mpc) 



140 
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Fig. 2. — Variation in the Hubble flow 6Hs = [Hs — H^^ / in spherical shells as a function of weighted mean shell distance: 
(a) CMB frame; (b) LG frame. In each case the fllled data points represent the flrst choice of shells in Tabled and the unfllled circles the 
alternative second choice of shells. We have omitted the flrst shell from the plots since 5H is so large in the CMB frame that it is off-scale 
- for the flrst shell: with a mean weighted distance of {r^} = 5.43/i-iMpc we have (^^CMB = ^-^^^ ^ ^-^^^ ^^LG ^ ^'^^^ ^ 0.007. 
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Fig. 3. — The Bayes factor \nB, where B is the Bayes factor for the ratio of probability that the LG frame Hubble flow is uniform in the 
region r > rs outside the shell with inner radius rs, divided by the equivalent probability for the CMB frame. 




Fig. 4. 



Angular variance of the Hubble flow in the full COMPOSITE sample: (a) CMB rest frame; (b) LG rest frame. 




Fig. 5. — Angular variance in the Hubble flow in the CMB rest frame for inner (r < To left panel) and outer (r > Vo right panel) spheres 
as To is varied over the values 12.5, 20, 40 and 60 /i-^Mpc. 




Fig. 6. — Angular variance in the Hubble flow in the LG rest frame for inner (r < To left panel) and outer (r > Vo right panel) spheres 
as To is varied over the values 12.5, 20, 40 and 60 /i-^Mpc. 
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TABLE 5 

Ratios C2/C1, C3/C1 of quadrupole/dipole and octupole/dipole for the multipoles of angular Hubble 

VARIANCE MAPS IN THE CMB, LG AND LS FRAMES, USING (fTTT) WITH NO IV VARIANCE WEIGHTING. In EACH CASE 
THE MULTIPOLE RATIOS ARE COMPUTED INSIDE (r < To) AND OUTSIDE (r > To) A BOUNDING SHELL. 



r < To { h~^M^c) 




< 12.5 


< 15 


< 20 


< 30 


< 40 


< 50 


< 60 


< 70 


< 80 


< 90 


< 100 


CMB C2/C1 




0.123 


0.061 


0.044 


0.098 


0.136 


0.191 


0.187 


0.167 


0.141 


0.134 


0.120 


CMB C3/C1 




0.010 


0.011 


0.007 


0.003 


0.005 


0.010 


0.009 


0.011 


0.011 


0.012 


0.012 


LG C2/C1 




0.653 


0.179 


0.123 


0.135 


0.116 


0.103 


0.104 


0.103 


0.092 


0.089 


0.085 


LG C3/C1 




0.067 


0.018 


0.008 


0.005 


0.008 


0.006 


0.007 


0.009 


0.011 


0.011 


0.011 


LS C2/C1 




0.861 


0.197 


0.133 


0.146 


0.124 


0.112 


0.113 


0.112 


0.101 


0.097 


0.093 


LS C3/C1 




0.068 


0.015 


0.006 


0.005 


0.007 


0.006 


0.007 


0.009 


0.011 


0.012 


0.011 


r > To ( /i~^Mpc) 


> 2 


> 12.5 


> 15 


> 20 


> 30 


> 40 


> 50 


> 60 


> 70 


> 80 


> 90 


> 100 


CMB C2/C1 


0.102 


0.096 


0.115 


0.124 


0.073 


0.038 


0.023 


0.041 


0.093 


0.093 


0.090 


0.327 


CMB C3/C1 


0.010 


0.009 


0.013 


0.015 


0.017 


0.009 


0.007 


0.011 


0.018 


0.078 


0.069 


0.076 


LG C2/C1 


0.072 


0.061 


0.064 


0.064 


0.053 


0.042 


0.032 


0.045 


0.068 


0.077 


0.066 


0.151 


LG C3/C1 


0.010 


0.009 


0.012 


0.014 


0.019 


0.013 


0.013 


0.014 


0.010 


0.048 


0.051 


0.016 


LS C2/C1 


0.079 


0.065 


0.068 


0.068 


0.054 


0.044 


0.033 


0.048 


0.074 


0.079 


0.070 


0.162 


LS C3/C1 


0.010 


0.009 


0.012 


0.014 


0.019 


0.013 


0.013 


0.014 


0.011 


0.053 


0.062 


0.020 



100 . 

CMB frame inside 
CMB frame outside 
LG frame inside 
LG frame outside 
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Fig. 7. — Ratios C2/C1, of quadrupole to dipole, on inner (r < Tq) and outer (r > Tq) spheres, as Tq is varied in the CMB and LG rest 
frames. 



Hubble flow variance and the cosmic rest frame 
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TABLE 6 

Least squares fit of dipole Hubble law ([Tq]) in radial shells in CMB and LG frames, for the 
TWO different choices of shells given in Table [T] For each shell, and each rest frame, we 

TABULATE Ng, Ts , fs (DEFINED IN TABLE [T]) ; THE BEST FIT DIPOLE HUBBLE CONSTANT, Hdg (UNITS 
hKMS~^ MpC~^, dipole slope (3s (units hKMS~^ MPC~^, THE GALACTIC LONGITUDE i^s AND 
LATITUDE bds OF THE DIPOLE APEX, AND THEIR RESPECTIVE STANDARD DEVIATIONS CF h ^ ^ , f^/3s, ^h^^ AND 
cr^^ ^ . We also tabulate the reduced (for u = Ns - 4) AND GOODNESS OF FIT PROBABILITY Qs IN 

EACH CASE. 



Shell s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Ns 


92 


505 


514 


731 


819 


562 


414 


304 


222 


280 


91 


Ts (/l~^Mpc) 


2.00 


12.50 


25.00 


37.50 


50.00 


62.50 


75.00 


87.50 


100.00 


112.50 


156.25 


fs (/i~^Mpc) 


5.43 


16.33 


30.18 


44.48 


55.12 


69.24 


81.06 


93.75 


105.04 


126.27 


182.59 


(^ds)cMB 


95.2 


87.2 


94.7 


97.5 


96.0 


97.9 


97.0 


99.7 


100.3 


98.4 


95.8 


(^^d JcMB 


1.3 


0.6 


0.7 


0.6 


0.6 


0.7 


0.9 


0.8 


1.0 


0.7 


1.9 


(^«)CMB 


95.8 


24.4 


14.0 


2.6 


4.9 


6.8 


3.8 


7.4 


5.8 


6.4 


0.9 


(^/3s)cMB 


1.5 


1.1 


1.2 


1.1 


1.1 


1.6 


1.7 


1.6 


1.8 


1.4 


2.9 


(^ds)cMB 


308.7 


300.1 


290.2 


309.6 


285.5 


279.3 


264.6 


274.7 


282.0 


309.0 


147.2 


V^^d JCMB 


1 Q 


o.U 


o.U 


oo.y 


±o.O 


19 7 


ZO.O 


1/17 


91 A 








-5.9 


17.3 


-5.8 


-42.0 


-3.0 


2.2 


12.9 


18.8 


11.7 


-2.8 


48.9 


(^^d JcMB 


0.9 


1.7 


4.5 


16.8 


10.7 


8.0 


20.5 


9.6 


15.6 


9.7 


195.9 


(XV^)CMB 


56.7 


6.1 


2.3 


1.5 


1.0 


0.8 


0.7 


0.7 


0.7 


1.0 


0.6 


(Qs)cMB 


0.000 


0.000 


0.000 


0.000 


0.430 


0.999 


1.000 


1.000 


0.998 


0.481 


0.999 


{Hds\Q 


79.1 


86.6 


94.6 


97.4 


95.9 


97.9 


97.0 


99.7 


100.3 


98.5 


95.9 


(^^^d Jlg 


1.5 


0.6 


0.7 


0.6 


0.6 


0.7 


0.8 


0.8 


1.0 


0.7 


1.9 


(^^)lg 


15.9 


20.0 


14.2 


14.9 


8.6 


4.9 


4.7 


1.5 


2.0 


4.7 


3.6 


(^/sJlg 


1.7 


1.1 


1.2 


0.8 


0.9 


1.0 


1.4 


1.3 


1.8 


1.4 


3.1 


(^ds)LG 


324.6 


52.7 


50.3 


87.2 


85.3 


86.3 


107.0 


265.4 


353.4 


348.1 


111.3 


(^^d Jlg 

(^rfs)LG 


12.7 


3.9 


8.3 


4.9 


10.1 


32.9 


27.9 


115.6 


217.1 


24.2 


76.0 


-30.5 


-34.8 


-55.4 


-36.3 


-44.2 


-58.7 


-42.5 


-51.0 


-76.2 


-38.1 


-17.8 


(^f'd Jlg 


6.4 


2.2 


5.5 


4.3 


8.0 


18.5 


20.7 


59.3 


57.3 


14.3 


36.5 


(xV-)lg 


7.8 


3.6 


1.6 


1.1 


0.8 


0.7 


0.6 


0.6 


0.7 


0.9 


0.6 


W^)lg 


0.000 


0.000 


0.000 


0.158 


1.000 


1.000 


1.000 


1.000 


1.000 


0.882 


1.000 



Shell s 


r 


2' 


3' 


4' 


5' 


6' 


7' 


8' 


9' 


10' 


11 


Ns 


321 


513 


553 


893 


681 


485 


343 


273 


164 


206 


91 


Ts (h~^Mpc) 


6.25 


18.75 


31.25 


43.75 


56.25 


68.75 


81.25 


93.75 


106.25 


118.75 


156.25 


fs (h-^Mpc) 


12.26 


23.46 


37.61 


49.11 


61.74 


73.92 


87.15 


99.12 


111.95 


131.49 


182.59 


(^ds)cMB 


87.5 


90.4 


93.8 


97.5 


96.1 


98.2 


98.5 


100.5 


97.2 


99.3 


95.8 


(^^d JcMB 


0.7 


0.6 


0.7 


0.6 


0.7 


0.7 


0.9 


0.9 


1.0 


0.9 


1.9 


(^«)CMB 


39.5 


16.2 


8.5 


3.0 


5.8 


5.2 


5.0 


8.9 


6.7 


5.4 


0.9 


(^/3s)cMB 


1.2 


1.1 


1.2 


1.0 


1.3 


1.6 


1.6 


1.7 


2.1 


1.5 


2.9 


(^ds)cMB 


294.9 


296.1 


307.3 


273.9 


276.2 


270.1 


289.3 


279.7 


314.8 


308.8 


147.2 


(^^d JcMB 


2.6 


4.6 


10.1 


20.2 


12.8 


17.8 


21.0 


12.0 


17.0 


20.7 


389.1 


(^ds)cMB 


17.9 


15.2 


-13.4 


-16.6 


14.9 


-0.1 


26.8 


7.0 


4.3 


2.1 


48.9 


(^^d JcMB 


1.2 


3.5 


7.0 


15.4 


11.4 


10.2 


16.9 


8.4 


13.9 


14.2 


195.9 


(XV^)CMB 


12.1 


3.6 


1.9 


1.2 


0.9 


0.7 


0.7 


0.7 


1.1 


0.7 


0.6 


(Qs)cMB 


0.000 


0.000 


0.000 


0.000 


0.910 


1.000 


1.000 


1.000 


0.132 


0.999 


0.999 


(^ds)LG 


88.2 


89.4 


94.2 


97.4 


96.2 


98.2 


98.4 


100.5 


97.3 


99.4 


95.9 


(^^d Jlg 


0.8 


0.7 


0.7 


0.5 


0.7 


0.7 


0.9 


0.9 


1.0 


0.9 


1.9 




20.9 


15.9 


15.1 


11.8 


5.4 


5.1 


2.8 


3.7 


5.0 


3.6 


3.6 


(^/3s)lG 


1.6 


1.1 


1.1 


0.8 


1.1 


1.0 


1.5 


1.6 


2.2 


1.6 


3.1 




42.9 


51.9 


64.0 


94.5 


94.1 


103.7 


70.4 


286.4 


0.0 


360.0 


111.3 


(^^d Jlg 


4.5 


4.8 


6.0 


6.3 


17.2 


33.9 


43.3 


35.4 


27.1 


33.9 


76.0 


(^rf«)LG 


-35.5 


-31.8 


-43.4 


-38.0 


-38.8 


-57.8 


-33.8 


-37.0 


-29.5 


-36.7 


-17.8 


(^f'd Jlg 


2.7 


3.7 


5.1 


5.1 


14.9 


17.1 


31.4 


23.1 


16.0 


24.4 


36.5 


(xV-)lg 


4.3 


2.7 


1.4 


0.9 


0.8 


0.6 


0.6 


0.7 


1.0 


0.6 


0.6 




0.000 


0.000 


0.000 


0.988 


1.000 


1.000 


1.000 


1.000 


0.404 


1.000 


1.000 
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TABLE 7 

Least squares fit of dipole Hubble law ([T9)) using data outside a radius 

r > To, AS To is varied, IN BOTH CMB AND LG FRAMES. WE TABULATE Tq, 
N(r > To), THE BEST FIT DIPOLE HUBBLE CONSTANT, Hd (UNITS /lKMS~^ MPC~^, 

DIPOLE SLOPE ^ (units hKMS~^ MPC""*^, THE GALACTIC LONGITUDE id AND 
LATITUDE bd OF THE DIPOLE APEX, AND THEIR RESPECTIVE STANDARD DEVIATIONS 

(^Hd^ ^/S' AND <ji^. We also tabulate the reduced (for u = N - 4) 

AND GOODNESS OF FIT PROBABILITY Q IN EACH CASE. 



CMB frame 



To 


N 


Hd 




/3 






^^d 


hd 


^^d 




Q 


15 


4358 


96.0 


0.2 


6.1 


0.4 


300.2 


4.3 


6.1 


3.1 


1.57 


0.000 


20 


4158 


97.1 


0.2 


5.8 


0.4 


289.6 


4.7 


1.0 
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Hubble flow variance and the cosmic rest frame 
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Fig. 8. — The slope (3 of the hnear dipole relation cz/r = + ^ cos cf), as given in Tabled is plotted by shell in the CMB and LG rest 
frames. The filled (unfilled) circles correspond to the unprimed (primed) shells. 
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Fig. 9. — All data in a sphere r > Tq is fit to the linear dipole relation (|19|) in the CMB and LG rest frames. We plot the slope, ^, of 
this relation as a function of the radius ro outside which the data is included. Open data points represent fits for the CMB rest frame and 
solid data points fits for the LG rest frame. 



32 



LG frame ro > 15.0 h'^ Mpc, AH : 18.4 h km 8"^ Mpc-^ (TV = 4359) 
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Fig. 10. — Angular variance of the LG rest frame Hubble flow given by in the r > 15/i~^Mpc range (solid lines), overlaid with 
angular uncertainties aa (colour map contours: (a) with no IV weights; (b) with IV weights (|18|) . The red dashed circle indicates a la 
region around the maxima/minima. Blue crosses indicate the residual CMB temperature dipole poles. 



Hubble flow variance and the cosmic rest frame 
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Residual CMB temperature dipole T(Sun-CMB) — T(Sun — LG) 
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FlG. 11. — Residual CMB temperature dipole in rest frame of Local Group (panel (a)) compared to LG Hubble flow variance dipole for 
r > 15 h~^Mpc with IV weightings (panel(b)). 
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Wiltshire, Smale, Mattsson & Wat kins 



TABLE 8 

Correlation coefficient between the residual CMB temperature dipole in the LG and LS rest frames, as compared to the 
Hubble flow variance in the r > 15/i~^Mpc sphere, for different choices of Gaussian smoothing angle a^. Unweighted (u) 

AND IV WEIGHTED (w) ANGULAR AVERAGING IS CONSIDERED IN EACH CASE. 



<T^ LG u LG w LS u LS w 

15° -0.8909 -0.9056 -0.8695 -0.8750 

20° -0.8945 -0.9197 -0.8774 -0.8965 

25° -0.8905 -0.9240 -0.8782 -0.9077 

30° -0.8847 -0.9237 -0.8769 -0.9133 

40° -0.8752 -0.9187 -0.8739 -0.9160 
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Fig. 12. — Pearson correlation coefficient for the correlation of the residual CMB temperature dipole in the LG and LS rest frames, as 
compared to the Hubble variance sky map in inner (r < To) and outer (r > ro) spheres, as ro is varied in the LG and LS rest frames. The 
correlation is computed both with and without the IV weighting. 



